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NATIONAL ADVISORY COMMITTED FOR AERONAUTICS 


TECHNICAL NOTE NO. 1018 


METHODS TOR DETERMINATION AND COMPUTATION OP PLOW PATTERNS 

OP A COMPRESSIBLE PLUID 



By Stefan Bergman 


SUMMARY 


A well-known method of generating stream functions of an 
incompressible fluid flow is that of taking the imaginary- 
part of an analytic function of a complex variable. In pre- 
vious publications of the author this method was generalized 
to the case of subsonic flows of a compressible fluid. Plow 
patterns, which until the present, have proved impossible to 
obtain by existing methods, were, however, obtained by this 
procedure; for example, flows around an obstacle the boundary 
of which is a closed curve, 1 as well as around nonsymmetric 
profiles. The procedure can be extended to the case of par r 
» tially supersonic flows. As this method for obtaining flow , 

patterns of compressible fluid from analytic functions of a,, 
compressible fluid requires rather lengthy computations, the 
« present paper is devoted to a detailed discussion of perform- 

ing these computations. The operations are divided into two 
groups: namely, those which need only be carried out once _ hTnd 
for all and then can be tabulated (or put on master cards) 
and those which have to be repeated In every individual case. 
A detailed description is given concerning the performance of 
necessary computations on punch card machines. This descrip- 
tion is illustrated by an example. 

In the appendixes some theoretical questions, which to a 
certain extent complete the results of NACA Technical Note 
No. 972, are considered. For instance, in appendix II, some 
questions which arise in connection with the determination of 
flow patterns around a nonsymmetric profile and the use of 
linear integral equations for constructing flow patterns are 
discussed. 


1 The method of Von Harman and Tsien yields a flow around 
a closed curve. However, this method assumes a linear pres- 
sure-specific volume relation, i.e., that p = A + (j/ft where 
A and o' are constant instead of the actual relation 
p = <Jpk (adiabatic case). 
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In appendix III those alterations are indicated which are 
necessary in order that the operator which has teen introduced 
for- subsonic flows may he transformed into an operator which 
generates stream functions of supersonic flows from two func- 
tions of one real variable. 


INTRODUCTION 


The mathematical theory of steady two-dimensional flows 
of an incompressible fluid is based essentially on the fact 
that a stream function of a flow of this kind can be obtained 
by taking the imaginary part of a conveniently chosen function 
of one complex variable s = 8 + 1 log v, where v is the 
speed (at the point) and 8 the angle which the velocity 
vector (at the point) forms with a fixed direction. 1 

The success of this method in dealing with problems of 
the theory of an incompressible fluid, seems to suggest the 
possibility of generalizing this approach to the case of a 
compressible fluid. An attempt, in this direction, has been 
made by the author in previous publications. To this end, in- 
stead of log v, there is introduced X(M), a function of 
the local Mach number M. 3 Further, instead of taking the 
Imaginary part of an arbitrary analytic function (i.e., apply- 
ing the operator Im(=Imaginary part of)) as in the case of 
an incompressible fluid, it is necessary to apply a general- 
ization of this procedure to obtain from f(8 + iX(K)) the 
desired stream function. 


x The introduction of functions of the variable s (in- 
stead of the customarily employed functions of x + ly, (x,y) 
being Cartesian coordinates in the plane) causes some diffi- 
culties of a mathematical nature; however, in contrast to the 
latter method, the former, more complicated approach (often 
called the hodograph method), admits of direct generalization 
to the case of a compressible fluid. 

s The function h(M) is real if M < 1, and purely im- 
aginary if M > 1. Thus, 8=8 + iMM) is a complex vari- 
able in the subsonic case and a real variable in the super- 
sonic case. Note that in the body of the text 
A(M) - 18 = — i ( 0 + iMM)) Is employed. 
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One of the advantages of this approach is that it mani- 
fests a far-reaching analogy 1 with .the case of an incompress- 
ible fluid, and is capable of yielding flow patterns which 
have not been obtained until the present - for example, flows 
around a closed profile, and so forth. 2 ThiB approach makes 
it possible to determine a flow pattern corresponding to any 
given function. In general, the actual construction of the 
flow leads to a considerable amount of computation; conse- 
quently; the use of. special computational devices such as the 
differential analyzer, punch cards, and so forth, would seem 
necessary as well as the preparation of certain tables which 
are independent of the specific flow pattern, and therefore 
need be prepared only once. 

The most efficient means of accomplishing this is not at 
all evident, and it is necessary to analyze the needed compu- 
tations from this point of view. The present report has been 
prepared in an effort to answer this question, especially as 
regards punch card .machines . 3 

This investigation, conducted at the Brown University, 
was sponsored by and conducted with the financial assistance 
of the Hational Advisory Committee for Aeronautics. 

The author was assisted by Mr. Leonard Greenstone, whom 
he would like to thank for his valuable aid. 


D1 SCRIPT I OH OS' METHOD 

In the author's previous report a new approach in the two- 
dimensional theory of a compressible fluid was developed. 


■••This analogy often serves as an indication of the proper 
method for obtaining results in the theory of compressible 
fluid which are similar to the case of an incompressible fluid. 

3 In the only case of a flow around a closed body which 
has. heretofore been considered, Von KArmAn and Tsien have as- 
sumed a linear pressure-specific volume relation, p = A + cr/p, 
A, cr being constant, instead of the actual relation p = crpk, 
k ® 1.4, which is used in papers of the author. 

The author would like to point out that other devices, 
in particular, the differential analyzer, are also of consid- 
erable importance for many of the above computations. (See 
reference 1 . ) 
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This method of attack is a general! eat i ojt of a proced- 
ure ordinarily employed in the theory of an incompressible 
fluids namely* the generating of stream functions of flows 
from analytic functions of a complex variable. 1 

Jcr the convenience of the reader the general idea of 
this method will be described in the following. The stream 
function \j/ of an incompressible fluid flow is a harmon.io 
function - that is* it satisfies Laplace's equation 

+ ili , o (i) 

Sx a By 8 

x,y being Cartesian coordinates in the plane of the flow. 
Conversely, every function which satisfies equation (l) may 
be Interpreted as a stream function of a suitable flow. Thus, 
if the Imaginary part of an analytio function f(s) of the 
coB’olex variable z = x + iy is taken, a stream function of 
a possible flow of an incompressible fluid is obtained. As 
eot 3-1 before, the method of generating stream functions in 
thin simple form cannot be extended to the case of a compress- 
ible fluid, since in the latter case the partial differential 
equation which \|/(x,y) satisfies is a very complicated non- 
linear one. This situation makes it necessary to use an al- 
ternate method, the so-called "hodograph method," - that is, 
to consider the stream function not as a function of s, 

but as a function of the velocity vector. 

If v lt v 3 , and (v,9) denote the Cartesian and polar 
coordinates, respectively, of the velocity vector v that 

is, if v a v x + iv a a ve 1 ® and if the stream function '\j/ 
is considered as a function of (v lf v s ) or of (log v,6), then 
\|/ is in each case a harmonic function of the given variables. 
That is, if \|/(x,y) is transformed by means of the substitu- 
tion 


x = x( v x , v s ) , 

y * y( v x .r 8 ) , 


5 ( x , y ) 

STvTTvsT 


4 o 


(2a) 


order to make possible this generalization, it is, 
however, necessary to consider the stream function in the so- 
called "hodograph" plane (i.e., in the plane the Cartesian 
coordinates of which are the components of the velocity vec- 
tor) instead of considering it in the "physical" plane (i.e,, 
in the plane of the flow). 
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or 


x = x(log v,8) , T ) 

+ 0 (3b) 

y = y(log V(0 ) i 3((log v),0) 

then the functions •*)/ obtained by the transformation (2a), 
(2b) satisfy the equation 

£±. + -lit = o (3a) 

B‘Vx a 3v s 2 

in the first case and 

— — +**Jl * o- ■ Ob) 

B (log v) s B0 a 

in the second case. (Note that these \|/* s are different 
functions of their respective arguments, although the nota- 
tion does not indicate this.) 

By writing 


\i/ = Im[g( v x - lv 3 )] 


(4a) 


or 


\J/ » Im[h(log v - i6)] (4b) 

if g and h are arbitrary functions of the complex vari- 
able Vj - iv 3 , and log v - 10, respectively, then the 

stream functions of possible flows of an incompressible fluid 
are obtained. 

Since the flow pattern in the physical plane is of pri- 
mary interest, it is necessary in this case, to carry out the 
transition to the physical plane; that is, to determine 
as a function of x,y. 

It is this second method which, though more complicated 
than the first, has the advantage of being capable of gener- 
alization to the case of a compressible fluid flow for which 
the equation of state, p * A + crp^, holds (A, cr, fc are 
constants, p the pressure, and p the density). In the 
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equation of state fsr an adiabatic process A = 0, however, 
this additional constant does not entail any theoretical 
difficult! e s . 


Aa has been indicated previously, the stream function 
of a compressible fluid flow, considered as a function of 
(x,y) - that is, in the physical plane - satisfies a nonlin- 
ear partial differential equation. If, however, is con- 

sidered in the hodograph or in the logarithmic plane (i.e., 
as a function of (?! ,r a ) and of (log v,6)» respectively), 
then \|/ satisfies, in each of these planes, a 1 1 near partial 
differential equation. 


In order to simplify this equation it 1 b expedient to 
introduce, instead of log v, a new variable A, 


8 106 


a l/3 a. l/2 l/h 

1 - (l - M ) ' /l + h ( 1 - M 3 ) 7 ' n 


1-1 + (1 


- M ) / l + h(l - H~? _ \ 

_ M a ) 1 ^ 3 n - h(i - k 3 ) 1 / 3 ' 


(5) 


where 


and 



here k 1 b the ratio of specific heats of the gas (k a 1.4 
for air), and a 0 the velocity of sound at a stagnation 

point. The equation which \J/ satisfies then assumes a par- 
ticularly simple form: namely, 1 

l 0 (\io s 1 falsi/ + SjJr) + H M 

4 V 3A a 3 0 a ' 3A 


0 

3z3z ' ' 3z 3z ' 


1 In the following, _instead of A and 0 the complex 
variables Z = A. - 10, Z = A + i 0 will frequently be used. 
The derivatives with respect to Z and Z have the follow- 
ing meaning . - - 


-9- » 1 ( 



a 1 / 

''_3_ 

i 

3 s 

„ 1/ 

/ 3 s 

, 3 a \ 

dZ 2 ' 

, 3A 1 

ae/ 

dz m 2 ' 

s. 3A 

ae/ 

3Z3"? 

4 

v3A a 

30 s / 
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where 


(k 4- l)M* 

8(1 - M S ) 3/S 


(7) 


In order to obtain a generalisation of the representation 
(4b), the author in the previous report derived the follow- 
ing re suit : 

From. the function N (see ( 7 )), certain other functions 
H(2X), ^ m ^(2X), n as 1,2,. . ., m = 1,2,. . . were deter- 

mined, and it is proved that the expression 



(where g(Z) is an arbitrary analytic function a, an arbi- 
trary non-negative constant) is a solution of (6). Thus, 
from an arbitrary analytic function it is possible to derive 
a function ^[X( v) , 8] , which represents the stream function 
of a possible (subsonic) flow of a compressible fluid. 6 Jor- 
mula (8) oan also be written in another form which Is suit- 
able for certain purposes: namely. 


1 It has been proved, subsequently, that it 
to interchange the summation and the passage to 
equation (8) to obtain 

z 4-, 


is possible 
the limit in 


(X , 0) * Im 


H(2X) [g(Z) + Y / • -J 

n=i 2 n ‘ 0 0 


c(£ 

n n 


•dtl 


(C 


However, as it is desired to make no reference to unpublished 
results, all the computations in this report are presented in 
such fashion that the use of (8) instead of (8a) entails no 
additional computation; equation (8) is almost always employed 
throughout the following. 

*l'ormula (8) can be considered as a direct generaliza- 
tion of (4b), since by choosing H = 1, for 

n and m, and Z = log v - IQ, (8) becomes (4b). 
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\|f(X,e) "» In| H( 2 'A) 


+1 


(^li) 


dt 


-i 


Jl - t' 


+ lim 
m-»-co , 


1+1 / 
S m (A.t) f (. 


Z(1 - t s > 


) 


dt 


t=-i 


yi - 


a - 1 


1 

J 


(9) 


E 


m (X,t) = 1 + y t Sn qi n) (2X - 2a) 


n=i 


where f(Z) is again an arbitrary analytic function of Z . 1 


(n) It; should be em P h asieed that the functions H(2A), 

^ ; (2X) are independent of the function g, and hence once 

computed (for a given value of k) may he employed in all 
other stream computat i one without change. 


Once 6) (corresponding to a given function g) has 

been computed, the transition to the physical plane - that is, 
the determination of the corresponding flow pattern in the 
physical plane - does not involve any;' thebr et i cal difficulty. 

Two problems immediately arise in connection with this 
method of attack. . ; 


I. How to determine, function, g, in (8), so as to ob- 

tain, in the phy si cal . plane , a flow, around a given obstacle 
or in a channel who se,. hpundgry curves ■ a^e, given . 



■ L "' — r* ** — tr r-e r- S ' j-fr i ’ — . 

Function f(Z) is connect ed ' wiili' ; g(Z)‘ by the follow- 

la nn : ■ t .. 


ing relation: 


f(Z) = - 


n v/a 

- Z sin d M-i? 2 "in h. gM 

n J d(Z sin* d) rr 
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II. Assume that g(Z) is known, 1 to develop a procedure 
which would permit the determination of the corresponding 
flow pattern in the physical plane with the minimum of compu- 
tation. Naturally, the flow patterns in which the aerodynam- 
icist is primarily interested are partially supersonic ones. 
Since the subsonic case serves as a basis for further devel- 
opments, as outlined in reference 3, 2 the author will limit 
himself in the present report primarily to this case. 3 

Although problem II does not entail any theoretical dif- 
ficulty, it does involve a very considerable amount of numer- 
ical computations for applications, as can be seen from the 
example described in reference 4, section 3, a fact which rep- 
resents a serious obstacle for the application of the method. 

Since the determination of various flow patterns is one 
of the purposes of the theory, the above-described situation 
suggests two possible modifications of the procedure for gen- 
erating flow patterns. ‘ 

1. The modification of the method so that a substantial 
part of the computation is independent of the particular 
choice .of g; thus these computations can be carried out and 
tabulated once and for all. 4 


x It may be remarked here that often a first approxima- 
tion to the desired flow pattern of a compressible fluid is 
obtained by substituting for g(Z) in (8; that analytic 
function the imaginary part G(log v,9) = Im g(log v - i9 ) 
of which gives the desired flow pattern in the physical plane 
for an incompressible fluid. The corrections which are neces- 
sary In obtaining a better approximation, as well as other 
methods of determining g(Z), will be discussed in future re- 
ports. (See. also reference 3.) 

s In sec'. 16 of reference 3 a procedure is described which 
makes it possible to generalize this method to the case, of 
partially supersonic flows. 

jj 

The author intends in a succeeding report to consider 
analogous questions for the case of a mixed flow in the light 
of methods described in sec. 17 of reference 3. 

4 

The need of tabulating various functions which appear 
in the theory of compressible fluids has been emphasized by 
some authors. (See, for example, Garrick and Kaplan (refer- 
ence 5), where the Chaplygin solutions have been tabulated.) 
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2. The rearrangement of the remaining computations 
(which must be Repeated in every particular case) in such, a 
form that they can he carried out with a minimum amount of 
labor using a punch card macshine, 1 

The main purpose of the present report is the develop- 
ment of a method of det er mining • the flow patterns according 
to requirements 1 and 2. 

In four additional notes certain problems considered in 
reference 3 are developed further; these are of a more theo- 
retical nature. 

In appendix II, the author shows that by employing re- 
sults obtained from a consideration of the singularities of 
the solutions of (8) and applying the theory of linear inte- 
gral equations, it is possible to determine a flow for a 
given hodograph. In certain cases, solutions of this kind 
can be considered as a first approximation to the solution 
of boundary value problems. 

In appendix III methods are given for the construction 
of purely supersonic flows, which methods employ various inte- 
gral operator representations. 

The derivation of the complex potential for a JoukowBki 
profile Is given in appendix IV, while appendix I is devoted 
to the question of determining the and- L^ n ). 


NOTATION . . 

The following list of notation is to serve the double 
purpose of being both an index of symbols used in the present 
report and a collection of some of the formulas, used in pre- 
vious reports, to which reference is made in the text; how- 
ever, no claim to completeness is made in this respect. 


has been emphasized by Kraft and Dibble, certain 
aspects of this theory may be successfully treated by use of 
the differential analyzer. (See reference 1.) 
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u = 
z 

Bu _ 
3 z 

1 (fis 

2 va x 

- i 

BuV 

By/' 

d 2 U 

= 1 ( 

3 a u + 

3 a u\ 

_ 1 

3 zB z 

4 V 

5 x 2 

By 2 / 

4 

a = 

a 0 - 

— ( k— 1 ) v 3 
2 

l/3 

a o 

speed 

of sound 

at a 


u_ = Su, = 1 /3u + t BuA 
z 3"z 2 '•Bx 3y/ 


speed of sound (reference 3,equa- 
tion (28)) 


J n 


coefficients in the series expansion of T in powers of 
x (formula (45)) 


coefficients in the series expansion of T 
of x (formula (46)) 


-l 


in powers 


e base of Naperian logarithms 
e xp ( x ) = e x 


fi»fa arbitrary twice continuously differentiable functions 
of their arguments 

g(Z) an analytic function of the complex variable Z 

g(-i)( 2 ) = ag(Z) = *g (o) (z) 

dZ dZ 


g (o) (z) = g( Z) 

g^ u ^(Z) nth iterated integral of g(Z) (formula (10)) 



k constant in the equation of state p = A + ap^: The 

ratio of specific heate at constant pressure to con- 
stant volume; k = 1.4 for air (reference 3, formula 
( 22 )) 


1(H) = ( — ~~*^ S ( 1— M 3 ( H ) ) 

^p(H)/ ' BH' 


(reference 3, formula (45)) 



+ cl 
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p pressure 

p 0 pressure at a stagnation point 

schlicht = univalent 

v magnitude of the velooity vector v; 
reduced speed v/a 0 

■ > ■ ■ ^ 

v velocity vector; that is, v * ve 1 8 

v lt v a Cartesian components of v; that is, 


occasionally, the 


V 3 Vj + iv a 


(x,y) Cartesian coordinates in the physical plane 




x s e 


note that in reference 3, X = a 


aX 


i-h\ 1 / h 


- = (lit) 


A constant in the equation of state p 3 A + ap^ (refer- 
ence 3, formula (32)) 

A * ~( k+1 )M 4 (M 3 -1 )~ 3 (formula (67)) 

8 


A n coefficients in the series expansion of T in powers of 
Xj. (formula (46)) 

A n,m coefficients in the series expansion of g(Z) as a 
function of £, of the fractional powers of £ 
(formula (17)) 

+ «*• 

A n , m ~ n iax( A n ^ m , 0 ) , A n,m ~ max ( ” A n , m » ® ) 

B n coefficients in the series expansion of T*" 1 in powers 
of Xi (formula (46)) 

®n,m coefficients in the series expansion of g(Z) as a 
function of £, of the integral powers of £ 

( formula (17)) 

,m = max( On , m, 0) J °n,m ■ aax (-0 n f a , 0) 
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E ,E /See appendix 111, bsc. 3); (reference 3, 

J “ theorem (53).) Note that in reference 


e 3 


B * = exp 


( 


Nd 


>/ 

— 03 


(!+£)) ■. 


differing from 


1 '-(i s £ + 


the usage her e . 
(k+l)M 4 


. -( 3k-l )M 4 - 4(3-2k)M S + 16 

64 ( 1-M S ) 3 L 

(formula (42)); (reference 3, formula (7l)) 


?jb( 2X) polynomial approximation of mth degree in X! to 


G operation a computation which, since it is independent of 

the flow, can he computed and tabulated once 
and for all 


H = exp (- 



J+i 


Sfd (X+£) = 




, i / a ( k- 1 ) 


(l-M s )‘ 

for the subsonic case (reference 3, formula (ill)); and 


l+I(k~l)M s J 

2 


E 



for the supersonic case; 


(formula (68)); 


in this sense H is used only in the series expansion 
of \j/ , as formula (8) 

(formula (52)); (reference 3, formula (42)); in 
this sense H is used as an independent vari- 
able. • 

Im' imaginary part of - 



= KS# + ^ 4 ) * * If 

4 \ 3A a 30 2 / 3A 


(formula (6)); (reference 3, 
formula' (46) ) 


L ^ n ^ ( 2A ) = - - - H(2A)<i^(2A) (formula (22)) 

2 n ni 
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t[ n) - l2nliH(2K)fiii n) (2N) 

m 2 n a 1 




M 


local Mach number ; M » v/a *» 


a o* ” |Oc-l) T * 


(formula (6)); (reference 3, formula (3l)) 


* 7 * 


N 


(fc+1) 


M 




b i J + i (formulae (49) and (8)); (reference 3, 

3 6 V ' formula (108)) 

functions, independent of the flow, which occur in the 
series expansion of \Jf (8ee formulas (49) and (8); 
reference 3, formula (84).) 

q( q ) computed employing F m instead of F (See 1 \ 

etc.) 

R*" 0 ^ ss H —■ (formula (23)); (reference 3, formula (114), ff) 
dv 

R^ a ^ functions, independent of the flew, which occur in the 

series expansion of \|/ (formula (23)); (reference 3, 
formula (114), ff) 

g(n) computed employing F instead of F 

m ni 


8 (l-M 3 ) 

A 


■j 7 g- (formula {?)); (reference 3, formula (47)) 


ts - 4 


F dA (formulas (49) and (8)) ; (reference 3, for- 


mula (107)) 


— 00 




Re real part of 

S operation a computation which must be repeated for each 

individual flow pattern to be computed 

8 # 1 - J (formula (44)); note in reference 3, formula (161), 
b is used for 1 - I . 


S 


(n) 


real part of that is, g^ n ^ =» + il^ n ^ 


T (n) 


imaginary part of g 


(n) 
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T = /l-M 2 

polynomial approximation of the mth degree in Xj, to 

W*(^,6; X^°^,0^ 0 ^); a fundamental solution of formula (6), 

possessing logarithmic singularity, 

at + iS^ 0 *. (See also 

reference 3, section 13.) 

x<p ) (T ,e) (See formula (35).) 

Yi p) (v,e) (See formula (35).) 
m 


a, ?p>t-i jO’Sp real and imaginary parts, respectively, of the 
: ’ " coefficients of z(p)in the power se: 


power series expan- 


sion of g(Z); that is, g(Z) 


-I 


(a 


ap+i 


+ ict 'ap )Z 


(p) 


(See formula (36).) 


p=o 


a a a max(a rL , 0) , a n = max(-c^ ,0) 


P(M) = 


tan 1 J M 3 — 1 — 1 tan” 1 ( h *fvi S - 1 X 
h V i . 


( formula ( 65 ) ) 


£ * Z + log 2 (formula (16)). In the appendixes £ = \ - i0 
T| = -9 + P(M) (formula (64)) 


8 


— > 


angle which v makes with the real axis 


6 ^ o ) i 6 ( 1 ) . . . . values of 8 at mesh points for a "lattice 

computation (See sec. 2.) 


X - 5 log 


i - J Hi? 


X- 


h 7l-M 


L i ■*" J i— H® i _ t^/i-M® 


5 1/h 


) 


(formula (5)); 
(reference 3, 
formula (48 ) ) 


X 0 corresponding to the maximum Mach number of a flow in 

X (o) i\( l) , . • . values of X at mesh points for a "lattice 

computation (See sec. 2.) 
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£ » 9 + £(M) (formula (64)) 


density; p = p 


1 - fc-1 v 8 
2 a n 


1 i/k-i 


(reference 3, for- 
mula (25)) 


P 

P. 


cp 

op 

* 


* 




* = 


modulus of £ ; that Is, £ = pe *'* 5 (formula ( 1 9 ) ff) 
density at a stagnation point 

constant in the equation of state: p = A + crp^ (refer- 

ence 3, formula (22)) 

potential function 

argument of £ ; that 1 b, £ = pe^ 0 ? (formula (19) ff) 

stream function 

»T+£ 

NdC£+£) 

/' 


exp 


NdC|+£)^ \{/, for the subsonic case (formula 

(4l)); (reference 3, formula 
-ext ( 69 ) ) ; and 

t+T) 


exp 



A( 


s)ds ] \l/ for the supersonic case (formula 

' ( 68 )) 


A Laplace operator A\j/(x,y) =* — — ^ - & » 4 f ■ ° & ■') 

Bx S Sy 3 ^ ded’z // 


X - 10 


Z a X + 18 

A m P JlW dH that is, M * J J(h) ( reference 3, formula 
./ (45)); (See sec. 5.) 

A„ * M 


H 


(n) 


oH 


X (See formula (33) ff.) 

(See sec. 5 of appendix 111.) 
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Remark : Observe that- quite, frequently functions will he con- 
sidered in different- planes although the notation may not, in 
general, indicate this. Thus, given f(xijr) , let 


x = x( x 1 , x 3 ) , 


y = y ( r 1 , x s ) , 


b(*L_e|) j 0 

B ( x , y ) 


and obtain fCxCx 1 ,! 3 ), yCx^-.x 3 )) = f 1 (x 1 ,x 3 ). The super- 

script will he omitted and only fCx 1 ,! 3 ) .written, since the 
meaning will he clear from the context. 


ANALYSIS 

1. An Outline of the Method to be Developed in- . 
the Present Report 


A method of determining the stream function (in the.... 
physical plane) corresponding to a given function g, the 
basis of which method is equation (8), was given in section 
3 of reference 4, together with a numerical example. An out- 
line of this method has been given. in the introduction. 

However, a good deal of numerical work is enta.iled by 
this approach, and the amount of labor involved increases con- 
siderably for a flow the maximum velocity of which is close- 
to the velocity of sound. For this reason, a modification of 
the method which would cut down the amount of .comput at i on is 
de'sirable. A description of the proposed modification fol- 
lows: 

The domain, D (in the ( X, S)-plane), in which the func- 
tion g(Z) is to be considered, can- be divided into two dis- 
tinct parts D x and D 3 , ' defined as the subdomains in which 
X 0 < X < 0 and X < X 0 , respectively. (See fig. 1.) The 
number X 0 is a preassigned number which can be altered to 
suit the case 1 although, in general, it will lie somewheVe- 


x The choice of X 0 ' will, depend upon the condition's in 
each case; for the most part, X 0 must be larger than the 
maximal X coordinate of the regular points' of g(Z) in V. 
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"between X = -0.4 and X = -0.1, corresponding to local 
Mach numbers M = 0.65 and M = 0.85, respectively (M as 
defined in equation (5)). 

In D x the argument X varies over values which are 

near aero, and, as a consequence the series (3) will converge 
very slowly, necessitating taking into account a great num- 
ber of terms in order to obtain a reasonable degree of accu- 
racy. On the other hand, g of equation (8) and therefore 
f is regular in Dj. and can be represented there by a 
aeries development . \ - •- 

In the domain D s the values of X are much smaller 

and therefore only a smaller number of terms of equation (8) 
need be taken into account. On the other hand, in D s the 

behavior of g may be considerably more complex; for exam- 
ple, g may have singularities and be many-valued. There- 
fore, it will be assumed that . g is given by its numerical 
values on a sufficiently fine lattice, or by a number of 
series, each of which converges in some subdomain of D s . 

In Dj the function g can be represented by a power 
series development, and since the operator (8) is linear it 
is hence possible to prepare tables once and for all, which 
will facilitate, to a very large extent, the determination 
of the flow pattern (in the physical plane). This will be 
explained in detail In section 4. 

In order to determine the flow in the domain correspond- 
ing to D s , the procedure of section 3 of reference 4 may be 

applied. Since the computations are rather extensive, It is 
expedient to employ mechanical devices. This requires a cer- 
tain modification of the above procedure, which modi f I cat 1 on _ 
will be described in section 2. Thus, two methods for deter- 
mining the flow corresponding to a given g(z) will be de- 
scribed in sections'2 and 3. Both methods employ punch card 
machines; in addition, the second method presupposes that the 
computor has certain tables available which are independent 
of the particular flow and hence can be computed once and 
for all. - - = — - - - •*' - 

Remark ; The divis_i_o.n of D_ into two subdomains Dj, and D-a 
is not neoessarily to' take place along the line X = X Q , 

It will often be more convenient to subdivide D as 

1 0r, if more convenient, by a polynomial approximation. 
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indicated in figure 2, so that the tables which have "been 
prepared may be used for the largest feasible part of the 
domain D. 

Remark : In order to emphasize the character of a computa*- 

tion which is being performed - that is, wheth.er it is one 
of that large class which need be computed and tabulated 
only once since they are independent of the particular flow, 
or whether the computation involved holds good only for an 
individual flow - to every description of a computation .will 
be added the characterization "(G operation)" or " ("S opera- 
tion)" according to its membership in the former or latter 
class of operations. 


2. Description of the First Method for the Construction of 
the Stream Function of a Compressible Fluid Flow by- 
Use of Punch Card Machines 


In this section the computation of a subsonic compress- 
ible flow by means of punch card machines will be described. 
This procedure is a modification of the method of section 3 
of reference 4. 1 


As indicated in that report, the procedure was divided 
into three separate stages. 


I. Computation of the integrals 

_ Z n-i 

g (a) (z) = 


J 1 g(Z n )dZ n dZ a _ x . . . dZ a (10) 


and the derivative ikil 


(o) 


dZ 


.(o) 


where g' ; (z) — g(z), Z = X - 1 9, is an analytic 

functi on 


II. Construction of the flow In the (X,0)-plane - that is, 
evaluation of expression (8) 


*In sec. 2 of reference 4, it was assumed that only an 
ordinary computing machine was to be used in performing the 
operations described there. 
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III. Transition from the logarithmic plane to the phys- 
ical plane 

Step I . - Three different methods of evaluation' of 

g (n )(Z), n-0, 1, 2, . . . and of (dg ( 0 > ( Z)/dZ) will be 

given in the following; two of_ these methods employ punch 
card machines; the third uses graphical means. 

The first method is to be applied if the real and imag- 
inary parts of g(Z) are given numerically on a sufficiently 
dense sat of points C ) of the lattice. 

The second method can be used when the function g^°^(Z) 
is given analytically and can be represented in the whole region 
concerned by' several series developments 1 * 3 around conveniently 
chosen points. 

The third method is much less exact; it can be used in 
order to check the results obtained by one of the above- 
described methods. 

Z 

g (n+i)( z ) x P g(n)( Zi ) dZi (11) 


may be written in the form 


(*,e) 

g ( n+i ) ( 2 ) ss P (S (n >dA + T( n >d6) 

' 0 o ( ^ > & ) 


+ i 


(T^dA - S^d0) 


( 12 ) 


where 


g 


(n) _ s (n) + iT (n) 


(See eouation (20) of reference 4.) The right-hand side of 
equation (12) may then be replaced by the approximating sum 
(13). 

1 Thase series developments are not necessarily power 

series since g can have singularities in D 3 , i^e., branch 
points, poles, etc. 

3 0bserve that cp^ n \ of sec. 3 of reference 4 are 

replaced by T n \ respectively. 
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^(n+i ) ( z ) = 

k=i 

+ (k-l)AX 


X Q + (k-1 ) AX , 6 Q jAX + i 


T 


(a) 


k=i 


• 6 0 ax + J T (n) x 0 ,e 0 + (k- 1 ) 


AG 


A8 


s 

-‘I s 

k=i 


k=i 

(a) 


X 0 ,8 q + (k-l)AG 


A8 (13) 


(See (3l) of reference -4.) ' The terms AX,A8 denote tix'e 
directed distances between the meshes of the lattice (see 
fig 3) ; that is, they are positive if the integration pro- 
ceeds in a positive direction, otherwise negative. 


As 


~(X, 0) 

•7(0, 0) 


indicated above, the (approximate) integration 
is to be carried at first from (0,o) to ( X, 0) along 


8 =0 (or if more convenient from (0,8) to ( X, 0) along 8 = 6 0 ) , 

and then from (X,0) to (X,8) along X = constant. 


A. (AH computations of A are (S operations).) The 


8 = 1,2,3, .. . can be 

computed on punch card machines by the following procedure. 
Every number 1 fs^°^[X 0 + (k-l)X,0 o ]j, k =1,2 is 

to be punched, say in columns 1 to 6, into a single card of 
a set Nj, . With every entry on this card an extra column c x 
(say, col. 7) is employed in which a number, say 1, is punched 
if S^°^[X Q + (k-l)AX,8 0 j is negative, and nothing is punched 
if the above number is positive. Then the cards are set for 
progressive totaling; (s^°^[X 0 + (k-l)AX,6 Q ]| will be added 
if nothing is punched in the column c x and subtracted if 1 

is punched in this column. The machine stops after each addi^ 
tion (or subtraction) punches the absolute "value of the pro- 
gressive total, in a new card s^, k = 1,2, . . ., n, say in 

l The symbol 11 | indicates that sign of 

be disregarded. 


sums f S (0) 


2^ S" J} j^X 0 + ( k-1 ) AX , 0 , 
k=i 


A X, 


has to 
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columns 1 to 6, and In an extra column c a , punches 1 if 
the total is negative and nothing if it is positive. 


Now the absolute value of AX is punched in an extra 
card M and, as before, 1 is entered in an extra column, c 3 
if AX is negative, and nothing if it is positive. Now, 
in a multiplying machine every number on the card s^ is 
multiplied with the number of M. In order to obtain the 
right sign, an extra column, e 4 , is provided in the new 
card. If the columns c s and c 3 , are both empty or both 

have 1 punched, then the machine will punch nothing into the 
column c 4 . If, however, in one of the columns c a (or c 3 ) 

the number 1 is punched and the other column, c 3 (or o g ), • 
is empty, the machine will punch 1 in column o 4 . 

The obtained results then have to be printed. In anal- 
ogous manner the remaining sums are to be evaluated. 


The obtained cards can then be used for evaluation of 
and T^ 1 ^, and so forth. 

Remarks : Clearly, the approximate summation can replace the 

integration, only if the integrand is uniformly continuous. 
Since, in general, the integrand has singularities, it is 
necessary to replace the approximate summation in the neigh- 
borhood of these points by the exact formula. This can be 
done, for instance, using series developments around the con- 
sidered singularity, (for details, see method B) or by other 
methods . 

The derivatives of first order, dg^°^/^2 * dg^°V^X, 
may be obtained by replacing differentials by finite differ- 
ences. (See method C.) 


B. A method, employing the example 


g(Z) 



2e “ ) 


zW s 


(i 



(14) 


considered in reference 4, which may be successfully applied 
when g( 2 ) is given by series developments, will now be 
described . 


Remark : If the function g(Z) it given in an analytic foxm, 
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then it ie always possible to represent it by finitely many 
series developments. 1 


In the case of the function given by the right-hand side 
of equation (14), the series given in reference 4, equation 
(25) may be used in order to represent g in the domain 
D a , A < -0.691. 


Another series development of equation ( 14 ) which is 
more suitable for the present purposes can be obtained in 
the following manner. 

The above function g(Z) possesses singularities 
(branch points of the second order) at the points 


Z = -log 2 + i ktt , K = 0, ±1, ±2, . . . (15) 


only. By classical results of the theory of 
can be expanded in series in powers of £l / 3 


funct ions 


g(Z) 


£ = Z + log 2 


(16) 


which series will converge for | £ | < 2rr and therefore will 
represent g in a large part of the domain, + D 3 , 

Wh i c h- Is of interest, 3 A formal computation yields 


= = i 


S = g 


U - sV /S + (1 - e £ r 1/s 


= i 


CO 

X. A °'" £ 


m- i/s 


^ = = i y 4 , 

dZ d £ - ’ 


m=o 


n=o 


n-i 


>( 17 ) 


S = i 



A 


,n + m- i/s 


n . m 


+ 1 ^ ^n , m ^ 


m 


n = o 


m=o 


1 A derivation of the analytic expression for the complex 
potential in the hodograph plane for a flow of an incompress- 
ible fluid around a Joukowski profile is given” inT appendix IT. 
By using this formula together with classical results of the 
theory of functions, the series developments for the above 
case can be derived. 


e Note that in example under consideration is deter- 
mined not merely in D s , but also in D x by the method de- 
scribed in the present section. 
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The values of A n>m and 0 n ,m are given in tables 5 and 6, 
respect ively . 

By writing 

g (n) c s ( n) + itU), n B _i, 0, 1 , 2, . . . (18) 


there is obtained 

CO 


n-i 


S (n) = - Y A nim p“ + ”-V=.in [( I1 + m-i) 1 ,j- Y °n,mP m •!» 


1 


me p 


(n) 


where 


m=o 

TO 

T - 

/ 

Bi=Q 


m=o 


n , m 


n+m- ijz 

p 77 cos 


m— l 


(n+m~i)cp + Y C nim p m cos m cp 


>(19) 


n=o 




£ » pe 1 ” 



The evaluation of 
chine proceeds as 

the S (n > and 
follows; 

T^ n ) on 

a punch card 

ma- 


The values of 

p^/ 3 , k as stl , 

±2, x3, , 

. .. p 1 /* - 

0.1, 

0-2, 

i • • t | of COB 

(|cp) and of 

Sin (|cp) 

k = ±1, ±2, 

• • • 

cp = 

0°, 30°, 60°, 

. , . , 330° can 

easily be 

computed ( see 


tables 7 and 8) and entered on three sets of punch cards A, 

B, C, respectively (G operation). By using set A, two new 
sets, D and E are then prepared (the following are all 
(S operations)). On every punch card of the set D the values 
of , m P n+m ~ 1 ' 3 and °f mP m ^ or a fi x ®d n and fixed p 

are entered, say A $ ~o p n ~ 1 / s are punched in columns 1 to 6, 

A + p n+1 / 3 in columns 7 to 12, A + p n+3 /® in columns 13 to 
18, and so forth. Here &n , m denotes A n>m if A n(ia is 
positive, and 0 if A n ,m is zero or negative; 0^ has an 
analogous meaning. In a similar manner A~ , ^p 11 " 1 " 111- ^ / 2 and 0 f 
c n,mP m ar ® entered on the cards of set 2. (Again A“ ( m = 0 
if ^n,m > and equals -A n ( m if H 0; the Bame holds 

f°r us i n S the sets 0 and D, 
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I 

i — r 


»+ n+m-i/s 
_p ' 

n , m“ 



is evaluated, and by using the sets 
computed 


n— i 

^ C n,mP m 3in mc P 
m=o 

0 and E there may be 



n , mf 


n+m- l / a 



n-i 

y °n,mP m 3in mc P 


( 21 ) 


m = o m = o _ 

By subtracting (20) from (21) is obtained. Similarly, 

can be determined. By interpolation, the values of 
S^ n ^(A,0) and T^ n ^(A,9) mey be determined at intermediate 
points. [Note that pe^P = (A - i 6) + log 2, which yields 
the relation between (p,cp) and (A, 6). Alternately, the ex- 
pressions (19) may be evaluated by adding on cards of the 
sets B and C an extra column, say, column 7, in which 
nothing is punched if the corresponding sine or cosine is 
positive and, say, 1 is punched if it is negative. In col- 
umns 1 to 6 the absolut e value of the sine or cosine is en- 
tered. 


Analogously, on cards of the set D an additional column 
is provided in which 1 or nothing is punched according to the 
sign of A n ^ m or C niin . The actual multiplication of the 

two factors proceeds similarly to that of method A. 

Since in the future it will be necessary to have values 
of S ( ^ and T^ n ) along lines A = constant, these val- 
ues for various values of 0 and for A = -0.02, -0.C15", 
-0.10, and so forth, were computed. 1 (See table 9.) 

C. The method described below is essentially the same 
as that described in method A; however, now, Instead of punch 
card methods, graphical means are employed. 

On millimeter paper the values of S^^A.B) and 

® ^ ° ^ ( A , 9 ) at first for some fixed values of A, say, for 

A = -0.02, -0.06, and so forth, and then for some fixed val- 
ues of 0, are drawn. (All operations of C are (S opera- 
tions).) (See figs. 8 to 15.) 

1 A portion of these values had already been computed 
(much less exactly) and presented in table II or reference 4, 
where the symbol T^ was used instead of T ( 11 > . 
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If S\(°\ T}^ Vj \ T g^° ^ are replaced by 

(AS^/AA), (AT (6) /A A )» (AS (6) / A8 >. (AT ^ 0 ^ /A0 ) , respectively, 

approximate values for \j/\ and tyg are obtained. (See 
table 10.) An integraph may^ be used to determine 

ys (6 ^(A,6)dA ( y^D (U ^(X,9)dX f / S ( o) (A , 6 )d9 .jAr ( (X , 9 )d0 , and 

forth (approx.), and so obtain, using equation (12), 

g (l >(z) * S (i )(X t 6) + iT^^A.S). Similarly, g (n )(Z), 
n = 2, IS, . . . can be computed. 

Ste p II .- The second stage of the method is trhen to ob- 
tain the values of the stream function and its derivatives 
in the (A,8)~plane - that is, to evaluate the expressions 1 

\|/(A,0) = L (o) (2A)T (o \a, 0) + L^ 1 ^(2A)T^ 1 ' > (A,e) + . . . 

+ L (n) (3A)T^ (A) + . . . 

L (cv) (2x) = H(2A), L (l) (2A) = 1 H ( 2 A)<} ( 1 > ( Z\) , . . . 

s 

L (n) (2x) = H(2A)Q (n) (3A) 

2 n n I 

\|/ v ( A, 0 ) « R (b) (2A)lm g 2 + R ( 1 ^ ( 2 A)T ( 0 } ( A , 0 ) + . . . 

+ R (n) (2A)T (n ' l) (A, 6) + . . . (23) 

l|/ e (A,0) » L (o ^(2A)Re g z + 1 ( 1 > ( 2A ) S < 0 >( A, 0 ) + . . . 

+ L (n ^ (2A)S^ n-1 ^(A, 0) + . . . (24) 

1 Since it is assumed in this case that the speed at 
every point of Da is considerably smaller than that of 
sound, the expression (S) is replaced here by 

\j/ (A, 0) = Im H( 2 A) e (Z) + V — | sll (2A)g < ‘ n ^(Z) 

- ■ L L 8 nl . J 

n=i 
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See equations (30), (31), and (33) of xeferenfee 4. Since 

the L ^ 6 ^ ( 2X ) , e => 0,1,2,. . . are independent of g, they 

caa be entered os. master cards once and for all, for differ- 

ent values of s and different values of X;. that is, they 
are (0 operations). E'er instance, on master card Ho. 1 in 
columns 1 to S, the value of L' & '(2X) for a fixed value X, 
say X 0 , is entered while nothing is punched in column 7 

if is positive; in columns 8 to 14 the absolute value 

of L^ 1 ^(2X) is entered, and in column 15 the number 1 is 
punched, if is negative, and so forth. . Similarly, -on 

master card No. 2 ths corresponding values ef L^ S )(2X( 1 )) 

are punched, and so forth. 


The remainder of step II consists • of ( 


From previous computation 
1,2,... for X = X l0? 


s the values 

0' 


s'- of (S operations) „ 

r of nfkx.eu k « 0, 

0' 1 ', 6 fa) , and so 


and 0. 

forth, for X s X^ ^ , 8 = , 8^ s \ and ,‘s* ‘fort-h, 

are obtained; both s-ets of cards, that is, the and 

T^ 8 ' arr then put into ‘the multiplier, which -then .yields -the 

values ©i 4 (22) for'the set of points (X^ 0 ^, 0^°)),. 

ej- 1 )},.- e< o) >. 

(s) 


( ) ’ ” " * ( * ) 

(X'“', 0 ° ), (X 3 , 0'"'),. . . And \Jz -7 -and‘ 4 ^g 'may- -be 
obtained, in similar- fashion . The-values -of ty, ty T , tyg ob- 
tained for the «as»~hnder~--c-on side rat i on are-given in tables' 
11 and 12. 




New, - the-value-S' of~' \J/(x(°) , 0^^), K - 0-,l.,2,. , ., 

( X 0 , 0 ) , \J/q(X 0 , ) are plotted on. graph paper' 

along the ab«eissa cf "which, the values of 0 are given. By 
using this diagram., the- values of 8 can be determined Tor 
which \J/(A^ 0 ^, 9) s*. -constant , say, 0 ± 0.1, =t0.2, and so 

forth. The values of \j/ T (-X^°\ 6) and of \j/g( x(<>) , 0) 

corresponding to \JKX^ 0 ^, 8) = constant may then be deter- 
mined. This procedure is then repeated for different values' 
of X.. 


See table 13 'and -figure 13. 

( 1C ) 

St en III .- To- every value X , K= 0,1,2,. . the. 

( KT ^ 

v * were determined for- which 


•values of 0 
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,.,YX ( n ) c ( n t ) \ 

^( A » o ) = T = constant 


( P5 ) 


as well as the corresponding values of (X , 6 ) and vj, (A,e). 

Tatles (or figures) of v 3 , 1 - M 3 , 1 <*A can 

pv 3 av 

prepared, which, since these quantities are functions of A 

alone, have to he computed only once, that is, they are (G 
operati ons). - • - — 

The image of a streamline (25) in the physical plane is 
given hy * 

- ■ - ^ 


x >= 


x(v) 


P 0 cos 6 j_ \Jr|(l - M s ) + V 3 v|f 

~n 


0 


pv 


3 

Z —= tv 


— 

y iB y(v) = - r'Po Bln & 1- - **> + v 3 j/ v a . 


>(26) 


dv 


Jo p vS ^ 

(See equation (19) 0 f reference 4.) J 

The integrals (26) will he appr o^imat ed _hy the sums 


x =* xl 


s = l 

J v i ) * ) Ax g ( t ) 

l—. 

S =1 

y = y ( Vj) = V Ay g ( T), l = 1,2,3, . , . 


>(27) 


Ax_ = 


P v 2 
s s 


p S= 0 

Ue (sT)c (1 - m s ) + v 


ay, 


Av c = 




s T » 




(st) 


+ v (st) £ 

( S t) 

cos 6 

S V j 

Vj/Q ( st) 

+ T*0 <ST>a 1 

a tr 1 

( s t) 

b in 6 

o V I 

~V ST) 


Av. 


Av. 


V — V 

S+3. s 


( 28 ) 
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The remainder of Bfcep III consists of (S operations). 
By -using tables for squares and the reciprocal, the values 


of 


*6 


( ST ) 


<r. 


( st) : 


and 


*8 


_1 _ 

( st ) 


are determined, together 


with the previously described tables for 
so forth. The quantity 

P° f . (st) 


pv- 


P* v 




(1 - K.) 


* .+ 


s s 


1 - M“ 


Av, 


and 


J \J/0 


( st) 


(29) 


is determined with the use of punch cards ; Equation (29) is 
then multiplied by cos 0*> ST ^ and sin 6' 8T ^ to yield the 
first and second terms of equations (38), respectively. 

Since the cosine and sine may vary in sign, an extra 
column must be provided with each term of the product as de- 
scribed previously. The cards are pv,; :n the multiplier 
which is set for pr 0|J e s R i.ve t the values (3?), 

which correspond to \2&) tnen resulting. ' ‘ 


3. Description of the Second Method for the Construction of 

a Compressible Fluid Flow 

As indicated in section 1, this method will often be 
applied if the local He ch number ie astray 1. ?cr this rea- 
son, in contrast to f\i' cci.oidereti o~.r? -nf section 2, it is 
now necessary to uss the exact formula 1 (8), that is, 

00 

' \ P r t ± ) ! 

12 t\i |g( * J-'-Iin ) 

1 L -i * o ? - n 

u 7U-VCC * n - Ct 

n =3. 

CO . 

= In L (o> C?,70g(r.} + lija V L m U 

m ->co *—> 

n = i 

z ln_i 

• • •/ s(£ a )a£ n - • • *£,. 

. 

L (0) (2A) = H ( 2 A) , L< n) (2)0 = H( 2 A)^ ( n) ( 2A) 

2 n i 

1 K method for det ermining the ihTTnT 

dix I . 



’(2X)g U) (Z)j. 


>J/(A,Q) «* lim \J/ m (A.6) = 


IrvH 

i 


th->oo 


is given in appen- 
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As indicated in reference 3, sections 9 and 15 (see alBO 
appendix I of the present paper), if X is considered only- 
in a range X < Xo < °. where Xo is & fixed negative num- 
ber, then a fixed m can he determined so that equation (6) 
can he replaced by 1 

5 ( *xx + W + *** “ 0 - (31) 


The solutions of-(3l) are given by 

n* i 

In the following it will be assumed that X 0 is a very- 
small number, say X 0 * -0.01 (i.e., that the flows with 

local Mach number = M Q >= 0,99 can be considered). Then m 
will be a very large but fixed number. 

Bemark ; In order to avoid confusion, all quantities which 
depend on m will have a subscript m; however, it is neces- 
sary to bear in mind that in this section m is a very large 
but fixed quantity, which remains unchanged in all considera- 
tions of this section. 

As Indicated in section 1, in tills method, certain tables 
can be employed which are independent of the flow and which, 
therefore, can be computed once and for all, and used in all 
subsequent computations. 


%,(X,6) = Im^L (o) (2X)g(Z) + £ 


*<Z)1 (32) 


1 Sirice N does not satisfy the hypothesis of theorem (8®, 
equation (6) was replaced by equation (31), where N m does 
satisfy the conditions of the above theorem and differs only 
slightly from N for values of X smaller than X 0 < 0. 

And A 0 can be taken as near zero as desired. 

In appendix I a method is given for determining N m for 
a given X Q with any prescribed degree of accuracy. 

Note that instead of -(\^xx + vJ/qq) + N m \{/^ =0 in appen- 
dix I, the equation. + ty* ) + * 0 is employed. 

This last equation is obtained from (3l) by means of the 
transformation (4l). 
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A. Description of Two Kinds of Tables 

1. If a sufficiently large number of the 1 . are 

computed, the functions X. m ( s p)(v,0) + ix m (®P + 1 ) (v, 0) which 

correspond to g(Z) = Z P ; that is, 
x m ( 2 p) (v,e)'+ iX< e P +l) (v,0) = H( 2 X)ZP 

OD 

+ V L^(2X)Z 

rfel 

where Z = A - i8 and A is given by (5), may be deter- 
mined . 


P+n L 


+ l). .. .(p + n) ( Gr operation) (33) 


Remark : In the case of 

A = log v, H ( 2 A) e 1, 
corresponding functions 


an incompressible fluid, where 
and L^ n (2A) «= 0, ■ n = 1,2,. . 
are 


the 


■x5 2 p) a Re (log v - i0) p_1 

^(ap + i) _ ( lc jg v _ iQ)P -1 


: (34) 


Analogously, as every function of 
every function x^ ' » p = 0.1,,. . 

and since for X < A 0 , N m practi 
one of these functions is a soluti 

2. To every function Xm' P ^ v > 
two real functions are determined: 


(34)'is a solution of (3) 
. . is a solution .of ,(3l), 
cal'ly equals N, every 
on 3 of (6). 

0), p — 0, 1,. • • 1 


9 


^See appendix I. 

Exactly speaking: An approximate solution of (6). It, 

however, does not differ essentially from the corresponding 
exact -solut i on of (6). 
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f P -s 

' (o,o) P 

(P )-) 


( 


-(l - M 3 ) cos 0 SX^ 5 

~a ae ’ 


sin 0 dXm 


dv + 


v (p) / c s 

Y m (v,e ) 


v a v J 

(v,e) 

• J (o,o) f 11 


cos 6 


axm 


(p) 


sin. 0 3X 


(p) 


m 


3 v 


39 r j 


de}; 


-(1 - M s ) sin 6 ax^p) 


>(35) 


cos 0 


ax£p> 


dv + 


sin 6 


Ip) 


av j 


cOCjgi 

3 v 


ae 


cos 8 3X 


ill] ae 


30 J 


} 




Remark : Since the above integrands are complete differentials, 
the values of the integrals are independent of the path of in- 
tegration (0 operation). 


Remark : In the case of an incompressible fluid, there is ob- 

tained for the corresponding functions X^p), Y^p) the ex- 
pressions: 


x ( 4p ) _ 

(P - 

I)" 1 

Pt'”- 1 ’ 

sin 

((p - 

r‘*> . 

(p - 

i r 1 

pr'*' 1 ' 

cos 

((p - 

( 3p + l ) 

(P - 

it 1 

(p-l) 

pv 

cos 

((P " 

( 3p + l ) 

SS 

(P “ 

i) _l 

(p-l) 

pr 

sin 

f<P - 



In the following it is assumed that the above-described 
funct ions, xW(r,e) and x£p)(v,0) y£p>(v, 0> are computed 

for a sufficiently large number of values of p and tabulated 
for a large number of values of (v,8). 


B. Determination of Plow Using the Above-Described Tables 

S_t§_p Ij Determination of streamlines ^hj(v, 6) = constant 

in the hcdograph plane: According to the assumption of sec- 

tion 1, the function g(Z) can be represented in the domain 
Dj, in which it will be considered in this section, in the 
form of a power series 
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g(Z) = ^ (a ap+1 + ia 2p )zP 

p=o 


(36) 


where ctp are real constant s , Jor s sufficiently large, 
the power aeries (36) can he replaced in by the poly- 

nomial 

I- <.*« * !.,>*’ < 37 ’ 

P~ o 


By substituting (37) into (30) and by observing (33), there 
is obtained for the stream function • corresponding to 

(37) 


^ m (v, 8) = Im 


y l/ P ^(2X) y (a„_+ia. 


) 


,n+p 


n=* o 
s 


)} 


2p Sp+1 . . , x 

p=o ( n+1 ) . . * ( n+p ) 

*« [I (4 sp> ^ sp+o y 

p= o 

S ^ p=2S+X 

- X (^ p xi Sp) ^ ap+1 4 3p+l) )= I (38) 


p = o 


p = 0 


Since, as a rule, it is necessary to determine the values of 
X m (v,0) at many points, it is convenient to use punch cards. 

3Tor every point (v,9) a master carql is prepared, and in 
this, in columns 1 to 6. the value of x|,o) at the consid- 
ered point (v,8) is entered and, in general, in columns 
6p+l to 6(p+l) the value of ( G- operations). 


Since a p can be positive and negative, (38) will be 
represented 1 by 

2S+1 ss+i 

( T • 9 ) = Z apX^ P ^(v,9) - Z a p X^ p) (v,8) (39) 

p=o p=o 

Each of the sums in the right-hand expression of (39) can be 
easily evaluated for a lj.r&e number of ooints using punch" 
card machines (S operations).. ?he curv e s v , 9 ) = constant 

can then be determined by interpolation. 

x a + = max(a,0), a - = max(-a, 0) 
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St e-g II.- Transition to the ■physical -plane. - To every 
point (v,0 ) of the h^dograph plane there corresponds a point 
(x,y) of the physical plane, which ie obtained by writing 


X =8 


>(v,e) 

*U.e) = / 

p l 


= y(v, e ) = 


(o,o) 
cos 0 

,( v J 0 ) 

\o,o) 


p 

aSa 

dv 

Po 


(l-M S )cos 0 o^ m 
v s 30 


sin 0 3\{/, 


sin 0 3\j/j 
v 90 


-^al dv 
3 v J 

- wj » ^ a p X^ P ^(v,Q) 

p = o 

°2S&I 




3 6 + 1 


>(40) 


f (l-M a )ein 0 3\|/ m cos 


.3 


30 


3v J 


'.in 0 5S5s + 

3v v 36 


ss + i 


I 


p = o 


a p Y m P)(v ’ 0) J 


See reference 3, equation (136); also equations (35) and (38) 
of this paper. 


bet \(/ m (v,6) ® e * constant be a streamline (in the 
hodograph plane). To every value of v on ^m(v,0) = c, 

there corresponds a value 1 of 8, say 0(v), which can be 
easily determined by interpolation or directly from the dia- 
gram for the tym(v,6).* constant. 


By interpolation (and.the ubb of th? tables described 
under II) the values of ^ ( v , 0 ( v ) ) , '(v,0(v)) are de- 

termined. Substituting these values Into (40) gives 
the coordinates ( x, y ) of the stream line \|/ m = c in the phys- 
ical plane (S operations). 


Remark ; Clearly, in order to apply this method, it is suffi- 
cient that the function g(Z) can be approximated in the 
domain under consideration by a polynomial 


s 



p*= i 


4* 


ia 


(e) • 
sp + i 


Z p 


On the other hand, by Runge's theorem, an analytic function 
can be approximated by a polynomial in every simply covered 

l 0r several values, say, 0 lt 0 g 0 n . 
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and simply connected domain. Because of this fact, the 
method described in this section may “be applied not only 
when is some region which lies inside of the circle of 

convergence of (36) hut for a much larger class of domains. 

CONCLUDING REMARKS 


A method of obtaining a subsonic flow pattern of a com- 
pressible fluid from a given analytic function g(z) is de- 
scribed in this report. The amount of time and labor needed 
for this method is reasonably small once certain tables have' 
been prepared. 2 These tables are completely independent of 
the flow, and consequently once prepared, the problem of de- 
termining the flow pattern may be regarded as solved, hot 
only from a theoretical but from a practical point of view as 
well. 


The present method yields only subsonic flow patterns, 3 
but by combining these with those described in section 17 of 
reference 3, it will then be possible to construct' mixe'd', - 
(i.e., partially supersonic) flow patterns, from a given 
function g(Z) . 


*The method described in this report , and references 8, 

2, and 3, is a generalisation of the determination of flow 
patterns of an Incompressible fluid from the complex potential 
g(£) = cp(log v,0) + I\j/(log v,8)£ = log v - 10, which poten- 
tial is given in the logarithmic plane. 

Assuming that the necessary auxiliary tables have teen 
prepared and that punch card machines, are available, the 
amount of labor needed in determining the pattern of a sub- 
sonic flow corresponding to a given function g(Z) will only 
slightly exceed that needed for determining the flow pattern 
of an Incompressible fluid from' a given g(lug v - i 6 ) . 

The author would like to emphasize that the tables of 
sec. 9 of reference 3, and those of the present report (the 
former are only an approximation to those of appendix I) serve 
merely to illustrate the procedure. The functions are com-^~ _ 
puted for comparatively few values of the arguments, and hence 
by using them it is possible to obtain only a rather inaccu«- 
rate picture of the flow pattern. 

jj 

Note that a similar method can be developed for purely 

supersonic flows. See appendix III. 
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A .method for determining various flow patterns is, of 
course, only the initial step in the study of compressible 
fluid flows, since the aer odynami ci st is, in the main, inter- 
ested in determining the influence of different factors such 
as the shape of the profile, the maximum Mach number, and so 
forth, on the flow pattern. 

By the choice of suitable functions for g, It Bhpuld 
be possible to obtain many cases of flows which are of con- 
siderable practical interest and value in studying various, 
phenomena in the theory of compressible fluids. 1 

Bemarlc : As has been emphasized in the introduction, it is ■ 

frequently of considerable importance to solve the "direct 1 .', 
problem determining the flow in the physical plane around a 
profile, which flow behaVes in a prescribed fashion at in- 
finity (i.e., far from the profile). Although in many in- 
stances it is possible* to determine the function g(Z) so • 
there is obtained a flow around a profile approximating the 
given one, it seems desirable to have a method of Bolving.the 
"direct 11 problem, and to determine when solution’s to this 
"direct" problem do or do not exist. The author hopes to re- 
turn to this question in a future report. 


Brown University, 

Prpvidenoe, H, I,, September 6, 1945, 


1 As has been indicated previously, the examples o.btained- 
which correspond to the Chaplygin solutions cannot, in 
general, yield the entire flow pattern (in the physical plane) 
around a closed profile. An exception to this has been the 
work of Itar man-T si e'n (references E> 7') but in order to ac- 
complish this they have substituted for the true adiabatic 
pressure - specific volume relation, a linear approximation 
to it. -- - 

S 0n ce a sufficiently large number of flows corresponding 
to various functions g have been "catalogued." 


NACA TN No, 1018 


3?', 


APPENDIX I 


THE DETERMINATION OP THE 



AND 



The operator (8) (see also secs. 9 to 11 of reference 3) 
was obtained in the following manner: As was proved in refer- 

ence 3, the function \}/(A,0) sati sf ie s' equati on (6) where.-iN 
is given by (7); A and M are connected by relation J 

Every solution \J/ of (6) can be written in the form 

ir « H(2M>lr*’ Ul) 


where H is given by (reference 3, (ill)) and \|/* 
the equation 


i f&hkl 

4 \ dA a 



P(2A)t* 


0 . 


sat i sf I e e 


(42) 


P(2A) 


1 t~ 5(l + k) 12k , + 2(3k - 7) 

8 L ( 1 _ M 3 ) 3 (1 - M 3 ). 3 1 - M 3 


+ 4(k + 3) - (3k - l)(l - 


M 3 ) 


(43) 


In order to determine P.C<?0 it is necessary to compute IT ■ 
as a function of X from (5)’ and ‘th'en- substitute into (43). 
The obtained function becomes infinite for A = 0 (i.e., 

for M » l), which causes certain difficulties. On the other 
hand, since only the subsonic case is considered here, and 
since a small modification of the function P(2X) practi- 
cally does not change (in the Bubsonic region) the solution »'f 
the equation 1 , it is expedient to approximate P(2A), in the 
range -<» X < A Q , |A Q I sufficiently small, by a function 

which remains finite at A = Q, for instance, by a p o 1 yr, omiai 
P m of the 2 m-th degree in a E ^. As was proved in reference 3, 
Tc (l-M ) ' can be developed in a series in e s N . i. fi nely; 

1 By u,sing the theory .of integral equations, it ie possi- 
ble to prove the following theorem. Let B be a given 
bounded domain in which; A A 0 V A 0 ’< ft and ‘in 'which T? m 

differs froM P by a suff icieno iy sm3.ll amount. To every 
solution of V(A, 6 ) ■ o'f 1 / 4 & 1 T + .P a 'l? + = O a solution ' 
t1/*('A, 6) of l/4Avj/* + .P>|/* = 0 can be so determined that 
]\|/* ( A, 6 ). - \{/ : ( A, 0) J A, 0 ) cB , and E i s a gi ven. small posi- 

tive number. 
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1 - T s x x + l( 2 k + l)^! + . . . 


(44) 


= 2 


( 


(k + l ) 1 / 3 - (k - I) 1 / 3 ./ k-i 


(k + 1 ) 


X / S + (k - 1 ) 1 /*' 


a aX 


This series converges for -e° < \ < o. Substituting k = 1.4 
into (44) yields 


and 


CO CD 

X A n y i ~ X a n x 


n=o 


n = 0 


CD CO 

-i " X Vi n - X b n xl1 


n = o 


n;=o 


(46) 


(46) 


Xj. = 0.239 8 


3 X 


x s o^\ 


The values of A n , a n , B n , b n are given in table 1 . 

Since (46) and (46) converge for -«» < X < 0; for -a> < X < \ , 
where X 0 <0 is a fixed quantity, it is possible to approx- 
imate (46) and (46) by polynomials 


m 

T m = X &n< 


sX n 


and 


n=i 

m 


(T 




n = 1 


V 


anX 


(4?) 


(48) 


By substituting these polynomials into (43) instead of 1 * 

. a g ( 1-M 3 )* 

and (1 - M ) , respectively, polynomials of approximation, 

3’m(2X) in e 3 ^ 1 , are obtained. Clearly, if a given degree 
of accuracy is required, m will increase as X 0 approaches 
0. By plotting I, l/T, T m and (l/T) n for a given m and 
comparing the corresponding values, the upper bound X Q of 
the values of X for which |P m (2A) - S’( 2 X)| is sufficiently 
small, may be determined. 
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TABLE 1 


n 

-A n 

— a n 

Bn 

Bn 

0 

-1 

-1 

1 

1 

1 

1 

.3392 

1 

.2392 

2 

1.9 

.1087 

2.9 

.1659 

3 

4.81 

.0658 

9.61 

.1315 

4 

13,939 

.0456 

33.869 

.1108 

5 

43.68 

.0342 

123.696 

.0968 

6 

144.02 

.0270 

462.39 

.0865 

7 

492.11 

.0220 


.0786 

8 


.0185 


.0724 

9 


.0158 


.0672 

10 


.0138 


.0629 


Bor instance, in the case under consideration where 
m = 10, the values of T and T m are given in table 3 and 
plotted in figure 16. As can he seen from figure 16, 

A < A 0 = -0.11 (i.e., M * 0.75), is practically 

equal to B(2A). If a good approximation is desired for 
bigger values of A, more coefficients 1 a n , b n must be com- 
puted. In order to check the obtained values of a n , as 
function of n, see figure 17. 

The coefficients Q,^ n ^(2A) of the operators which yield 

solutions of the equation \}/g ^ + B m v{f =■ 0 can be obtained in 

the same way as derived in reference 3, from which reference 
the results are obtained 

x It may be remarked that other methods of obtaining ap- 
proximating polynomials for B exist. These will not be 
investigated in the present report, despite the fact that 
they merit considerable attention. 
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Table 2 

The values of T, M, and Tio 


-2 A 

T 

10 

T 

nJi-T 2 = M 

0.0160 

0.43644 

0.300 

0.954 

0.0195 

0.44208 

0.320 

0.947 

0.0230 

0.44758 

0.336 

0.942 

0.0265 

0.45300 

0.350 

0.937 

0.03CO 

0.45835 

0.365 

0.931 

0.0335 

0.46360 

0.380 

0.925 

0.0370 

0.46877 

0.390 

0.921 

0.0405 

0.47385 

0.401 

0.916 

0.0440 

0.47885 

0.412 

0.911 

0.0475 

0.48377 

0.421 

0.907 

0.0510 

0.48861 

0.430 

0.903 

0.0545 

0.49338 

0.439 

0.898 

0.0580 

0.49807 

0.448 

0.894 

0.0615 

0.50268 

0.455 

0.890 

0.0650 

0.50723 

0.463 

0.886 

0.0685 

0.51170 

0.470 

0.883 

0.0720 

0.51610 

0.477 

0,879 

0.0755 

0.52033 

0.484 

0.875 

0.0790 

0.52471 

0.491 

0.871 

0.0825 

0.52892 

0.497 

0.868 

0.0860 

0.53306 

0.502 

0.865 

0.0895 

0.53714 

0.507 

0.862 

0.0930 

0.49225 

0.512 

0.859 

0.0965' 

0.49790 

0.520 

0.854 

0.1000 

0.54902 

0.525 

0.851 

0.1035 

0.55287 

0.530 

0.848 

0.1070 

0.53639 

0.535 

0.845 

0.1105 

0.56040 

0.540 

0.842 

0.1140 

0.56408 

0.545 

0.838 

0.1175 

0.56771 

0.550 

0.835 

0.1210 

0.57129 

0.554 

0.832 

0.1245 

0.57481 

0.559 

0.829 

0 . 128C 

0.57829 

0.563 

0.826 

0.1315 

0.58172 

0.567 

0.823 

0.1350 

0.58511 

0.571 

0.821 

0.1385 

0.58844 

0.575 

0.818 

0.1420 

0.59173 

0.579 

0.815 

0.1455 

0.59498 

0.583 

0.812 

0,1490 

0.59818 

0.587 

0.810 

0.1525 

0.60134 

0.591 

0.806 

0.1560 

0.60446 

0.594 

0.804 

0.1595 

0.60754 

0.598 

0.801 

0.1630 

0.61054 

0.601 

0.799 

0.1665 

0.61357 

0.604 

0.797 
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ft 
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m 
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APPENDIX II 


THE EQUATION (IN THE GANONIOAL POEM) POE TEE POTENTIAL FUNCTION 
AN APPLICATION OP INTEGRAL EQUATIONS TO THE 
THEORY OP COMPRESSIBLE FLUIDS 


1. In section 6 of reference 8 and section 7 ef refer- ■ 
ence 3 the equation (in canonical form 1 ) for the stream func- 
tion has "been derived. See equation (6.6) of reference 8 or 
(46) of reference 3, 

There are instances, however, where it is more conven- 
ient to operate with the potential function $ rather than 
with the stream function \}/. 

In this section the canonical form of the equation for 
cf) will he derived. 

1 By introducing suitable new variables | = £(x,y), 

T) = t|(x,y), every equation L(\l/) == aty xx +. J3b^ xjr + c^yy + d\J/ x 

+ = 0 of elliptic type can be reduced to the form 

^11 + ^T)T| + + + 0^ = 0, so-called " canonical form of 

equation L . " (See reference < 3 .) 

In the case considered in the present section i s H, 
y = 6 , and i = X, T) = Q . 
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Functions 4> and satisfy the system of equations 

( B4>/ 30 ) » (3\|//3H), l(H)(a^/de) * -(S45/3H) (50) 

De’quation (6. 21) of reference 8 and equation (30) of refer- 
ence 3] 

where 

(dH(v)/dv) * p / v , 1(H) * (1 - M 3 )/p® (51) 

[equations (6.l), (8.18) of reference g.;. equations (43), 

(43) of reference 3], 

If, now, the new variable X, given by 

1 A 

(dA/dH') = (l - M 3 ) /p, that is, (dX/dv) = ( 1 - M 3 ) /v (53) 

[equation (6,4) of reference B and equation (48) of reference 
3], is Introduced (50) becomes 

« 0 » P -1 (1 - M 3 )H x , p** l (l - M 3 )^v|/ e = -(J> x (53) 

Differentiating the first equation (53) with respect to Q 
and the second with respect to X yields 

4)qq ** p _1 (l - M *)^4/X0 (54) 

P _1 (l - M )^X8 * ^d(p _1 (l - M 2 )^)/dX ;\|/q * -<j>xx 

Replacing the first term of the second equation of (54) by 
4>00 and \|/ Q by -p(l - M S )“^4) X (see (53)) yields 

4 >ee + ^xx - p(i - M s ) _ ^fd(p _1 (1 - M s )^)/dx~|<!>x * o (55) 

U -J 

Now, by the second relation of (53) 

pd - M S r* ["d(p“ 1 (l - M a )^)/dX 

L 

= p 3 (l - M S )~ 1 (v/p)['d(p' rl (l - M 3 )^)/dvl - 4N 

L - 

= - <* t (1 _ 


(56) 
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[equation (6.6) and errata of reference 8, or equation (47) 
of reference 3], Thus, the" equation for 0 "becomes 

(57) 


4>XX + *ee - 4N<t^ = 4 !"<$»££ - BT(<*>J + <fc>£ ) 


t = X(v) - 16 


2. A case in which it is more advantageous to ' consider 
4> .rather than \J / is the following: : 

, ... In section 7 of reference 8 and in section 13 of refer- 
.ence 3 singularities of functions satisfying equation (6) 
were considered. As was indicated there, a, flow with a. 
"vortex-like" 1 singularity at , (X°, 0°) .is obtained if' for 
the . stream function, the so-called fundamental solution 

tf*(X,e ; , 5 ( 0 )) ' 

= A(Lf; £<*> J(o)) log \i - t C l + (58) 


[equation (7'.1) of reference 8; equation (119) of reference" 3] 

£ ® X - 10, £"=X+i0‘ . - - - r 

i s taken . 

A 8 was explained in section 14 of reference 3, it is Im- 
portant (in connection with the transition to thephysical 
plane) to have (working in the X, Biplane) singularities the 
derivatives of which with respect to X 'and- to -0 are 1 1 

single-valued functions of X and 0. - . ■ 

The point ■ £ 10 ' corresponds to the pdint z - 03 of the. 

physical plane, and if for the potential function, 4>, a- 
fundamental solution 1 ' . ' , w 7~ 

A*(C,Ti £ (o) ,£ r °l) log U - £ ? | +B»(£,r : £ (o) ,? r ° y ) (59) 

of (57) is taken, a flow with a " Bouroe-like" singularity is 
obtained, (l&xpre ssi on(59 ) and, therefore , it s derivatives are 
single-valued .fuhcti ons ' of X and 6.) 



The names "vortex-like” and "source-like" are used be- 
cause in the case of an incompressible fluid (and 'in the phyg - 
i cal plane ). . in the case of gi -vortex .the stream, function is 
given by m log is — z° I , and in the case of a source the po- 
tential function is given by 4> = m log Iz - z°f , m being a 
real constant. (See reference 10, pp. 19 p ?nd 320.) 
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A single-valued, solution of (57), which is infinite of 
the first order at £ = £(°) may be obtained by taking the' 
derivative with respect to Q of (59). 

3. A problem of considerable interest is that of deter- 
mining a flow of a compressible fluid around a given profile, 
or at least around a profile the shape of which approximates 
the given profile. Since, in many instances, by reasoning 
from the incompressible case, the approximate image in the 
hodograph or (A,0)-plane is known 1 , it is possible to con- 
sider, Instead of the above problem, the question of deter- 
mining a flow for a given hodograph . and the behavior at the 
point of the hodograph corresponding to z =» co 8 i b prescribed. 
Clearly, instead of the image in the hodograph plane the image 
in the (X,0)-plane may be used. ' If the results of section 7 
of reference 8, section 13 of reference 3, and those of sec- 
tion 2 of this appendix are employed, it is possible to deter- 
mine a function ^(^,0) satisfying (6), which possesses 

the required behavior at g = oo. Naturally, \|/i ( A , 0 ) for 

the point 3 (\»* Q^) must have, a singularity which satisfies 
the conditions indicated in section 14 of reference 3, in or- 
der that the flow in the physical plane will be a flow around 
a closed, curve. (See, in particular, equation (145;) of refer- 
ence 3.) Function v[^i(A,0) is as yet, not the required 
stream function, since it does not assume constant values on 
the boundary of the domain. In order to determine this func- 
tion, it is necessary to find a solution ^ S (A,0) of (6) 
which is regular in the domain H x , and which assumes, 
on the boundary h x ■ of H^, the values 4 

r The image in the logarithmic plane of an incompressible 
fluid flow around a profile P is often used as a first ap- 
proximation of the image In the ^8-plane ■ of the flow of a 
compressible fluid around a profile similar to P. See figs. 

4, 5, and 6, where the boundaries (and some streamlines) of a 
flow around a Joukowski profile in the physical, hodograph 
and (paBudo-) logarithmic plane, respectively, are given. 

2 The coordinate z refers to the physical plane. 

3 Th.e point (^,9^) corresponds to the point % * oo of 
the physical plane. ... 

4 

Since the domain Hi extends to infinity and, in gen- 
eral, is. multiply covered, it is necessary to alter somewhat 
the method of attack to be described, by mapping Hi conform- 
ally on a finite and Schlicht domain. 

For the sake of brevity this stqp will be omitted in the' 
following . ' " * ' ■ 
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WV = "WV ' 60 > 

“being an arbitrary point of h x . 

Function ta(^,8) can be determined using the theory of 
integral equations. (See footnote 8, p. 281 of reference 2.) 
Indeed, let t 3 (X,0) be that harmonic function which assumes 
the prescribed values on hj. , then 

^ 

satisfies the equation 

at 4 at. 

At, + 42T — ± = - 4N — ^ (61) 

O A O A 

and vanishes on the boundary h t . 

By employing classical results t 4 can be obtained as 
the solution of the integral equation: 


t 4 ( t.e) = *4 ( \.e x > 4 \ de x + 


t B = -3w JJ 4H 0- dX x dBj, 


x 

where G = G(X,0; X lt Gi) is Green’s function (of Laplace’s 
equation) with respect to the domain Ex. ' ~ 
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APPENDIX III 


A METHOD FOB. DETERMINATION OF STREAM FUNCTIONS OF 
PURELY SUPERSONIC FLOWS 


1. As indicated, in reference 8» section 10 and reference 
3, section 16, the approach developed in these papers makes 
it possible to construct mixed (i.e,, partially subsonic and 
partially supersonic) flows by use of the following procedure: 

In preceding papers two methods have been described 
(one given by Chaplygin, the other by the author 1 ), which 
yield certain types of particular solutions the stream 

function of a compressible fluid flow. (See sec. 8 (8.3), 
(8.6), (8. 22) of reference 8 and sec. 2 of reference 2.) The 

represent stream functions of flows, which, in general, 
include subsonic and supersonic fegions. 

As was pointed out in detail in reference 2, section 3 
and in the introduction of reference 3, the flow patterns 
generated by the \J/ V mentioned above or a linear combination 

of them are of rather special character. In partic- 
ular, the flow patterns with stream function cannot 

(in general) represent an entire flow around a closed body. 

Frequently, In the theory of analytic functions of a com- 
plex vgiriable in a similar situation (i.e., when one expres- 
sion of a certain kind - e.g., power series - does not repre- 
sent the function, say, f, in the entire domain B In which 
the function has to be considered), the procedure employed is 
to decompose into smaller regions, say, into B%, K a 1, 

2 ^ Bjj. = (see fig. 7) such that it is pobbI- 

ble to find in every region Bj^ , another analytic expression, 
Bay, f £ , which represents f in that region. G-enorali sing 

1 Bers and Gelbart, in reference 11, obtained the same 
solutions independent ly of the author. They denote the func- 
tions <py + I \|/ v as S-monogenic functions. Here qp v , is 

the potential funotion which corresponds to the stream func- 
tion vj.y. 
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this method of representation of a function of a complex vari- 
able, the author described in section 10 of reference 8 and 
in section 1 7 of reference 3 a method for representing the 
stream function and in a similar manner; that is, in decompos- 
ing the domain B into parts B K , and representing \|/ in 
every B^ by another analvt i cal expression. 

In order to apply this method 1 a representation for a 
purely supersonic flow is frequently required. 

A method for generating purely supersonic flows, com- 
pletely analogous to that developed for the subsonic case, ■ 
will be given in this appendix. 

2. The equation 

S(\|/) mf'BjLf (1 - M S ) = 0 (63) 

V p / 38 3 E s 

(equations (43) and (6.2) of references 3 and B, respectively) 
serves once more as the starting point for the following con- 
siderations. 

In order to write the equation for \j/ in the "canonical 
form" 2 * * , it is- necessary to introduce new variables t.T]. 

1 = 8+ p(JC), T1 = -8 + 300 (64) 


where 



A It may be noted that purely supersonic flow patterns can 
occur upon considering flows In channels or around a body with 
a cusp, in which case the flow has no stagnation point. 

2 Since in the supersonic case M > 1, equation (63) is 

of hyperbolic type. By introducing suitable variables ^ ,Tj . 

every equation of hyperbolic type can be transformed’ into l;he 

so-called canonical form + + 0 . 1 ■ ' 
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Equation (63) then becomes 

+ A(\J/| + xj/^) = 0 

where 


A = i(k + 1)M 4 (M 2 - l)“ 3/a 
8 


The function 


= H\|/ , H = exp 
satisfies the equation 


,( I -nn) 


A ( S ) d 8 


xl/t - F\j/* = 0, F = A s + (dA/ds), a = t + T] 


( 66 ) 


(67) 

>(. 68 ) 

(69) 


3. By use of considerations similar to those developed 
in references 8 and 3, the following theorems can he derived: 


The or em I . - Suppose that S' m is a function which pos- 
sesses a continuous first derivative. Let Bj;(S t T),t) and 
®&(5|T),1:) he solutions of 


(1 - t a ) 


d a n; 

6T| dt 


as I 

an 


3t i 




d a z* x 


a|an 


_ F S 
m i 


= 0 


(70) 


and 

(1 - t ) 
respectively. 


£lk 

a|at 


1 ® 2 

- — + 2tT| 

t a i 


2 * 

• 5 Ba 

.aean 




(71) 


Let B], and $ s possess continuous second derivatives, 
and let ( aE*/d| ) /T] t and ( SEg/ dT])/ he finite for t = 0^. 
Then 


u( |.n) 




[s*( i.n, t)f x (ild - t s ) ) 

+ E*(j,Ti,Of 3 J.T](1 - t a ))J (1 



(1 c 


(72) 


where fj^, K = 1, 2 are two arbitrary, twice continuously 

differentiable functions of their respective arguments, is a 
solution of the equation 
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.. ip U = 0 (73) 

agan m 

The proof of this theorem is given in reference 12 > section 2. 

Theorem II .- Let r m (£) possess derivatives of all or- 
ders in- the interval P.q Is ^ = (3 1 » 0 < P 0 < < ®* Ifa 

constant a exists such that the inequalities 


a Kg m ^ c(K + 3) 
dp K + 3 ’ 


0 1 1 1 2 j . 


e 0 ^ p ^ 


obtain, then there exist solutions 1 1 ( | , T| , t ) and B s (|,Tl,t) 

of (70) and (71), respectively, satisfying the conditions of 
theorem I . 

By substituting the functions Eg, K = 1,2 into (72) 
for the E k , there is obtained a representation for solutions 
of equation 

- V* = 0 (?5) 

in terms of two arbitrary, twice differentiable functions. 


' 4. There' exist various other integral representations of 

solutions of (69) in terms of two arbitrary functions of one 
variable. One such representation, differing from that given 
in the preceding section, will be discussed here. 

Let & (| , T) ; | * , T|* ) denote the ftiomann functions of equa- 
tions (69). (see reference .9., p. 22) - that is, a function 
of the four real variables \ ,Tj , £* , T1 ^ , which satisfies equa- 
tion (69) for every fixed (|*,T|*), and which further has the 
properties that. 

ft(l,T] ; I ,r\ ) * i 

\ ■ (76) 

, . * . * # . . I 

ft(| ,T);£ ,r\ ) = 1 
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This function (for (69)) may he . repre sented in the form 


% ^ 


B . *0 ; I ■**“, Tl* ) = 1 - 



+ / / J’(| 1> T1 1 ) 

(See reference 3, sec. 7.) 


S’(i x .Tl 1 )d| 1 dT) x 


2 , (e 3 ,n s )de 1 dT], &i 3 dn s 



.(77) 


The classical theory of partial differential equations of 
hyperbolic type yields the following results; 


Let fj[(K * 1« S) he any two arbitrary differentiable 


functions of one real variable, and if u satisfies the dif- 
ferential equations ug^ + Fu = 0^ then 


u (| , T| ) =» u(0, 0)R(g ,T) 50 , 0 ) + / R(S t 11;g* l 0)f 1 (!*)d6* 


o jn 

+ / R(£ ,T1;0,Tl*)f a (T1*)dTV» (78) 


(R is, of course, the Riemann function associated with the 
differential equations satisfied by u.) 


Remark ; A representation of the form (78) is a? 
the subsonic region. 


j. i :i £or 
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Indeed, suppose that | is replaced by £ = X - i8, T1 by 
X = X + i$, t* by 1* , Tj* by f* , in (77). this trans- 
formed expression will not differ essentially from the func- 
tion x( S.D introduced in (7.4) of reference 8. 

If the complex variable £,"][" instead of X,6. is used, the 
equation for the function , in the subsonic case, ass'ufieV 
the form" 


(equation (86) of reference 3) 


a n8w representation for in terms, of two arbitrary ana- 

lytic functions h lf h s of one complex variable may then be 
obtained: £ 

VU.T) - **(o,o)h(£,T;0.o) + f. a(S.Tsl f .o)h x (£*)dS* 

79) 


(R. is the Riemann function of the differential equation 
for V 4 .) - - — 

5. It is of considerable interest to show that both (8) and 
(72) are different forms of the same operator, the former ob- 
taining in the same subsonic case while the latter holds in 
the supersonic case. In order to derive this conclusion. It 
Is necessary to develop further the method of attack initiated 
in sections 6 and 8 of reference 3. The following result is 
a slight generalization of theorem (53) of reference 3. 



Ra,£ ; o,r)h 2 (£*)d^ ( 


Let E 1 be a solution of 

Att 3 ( — 
H \ Al 


o^u.e.n s [A - 


l_t (A + 10) 


3E 3® x Ahh 3 i 

Afl SE + 1 3 0 + 2 A h S 


) 


. fX 

36 s 


3 3 2 

VP 


f- + = 0 1 (80) 


Hote that A is a function of H alone and that 


A H = |A. A hh = Indeed, A H 3 = l(H) = p 0 S p- 3 (l _ M 3 ) 

and A(fe) = X (M) . (See sec. 8 of reference 3.) 
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which possesses the property that 


( i a h a vi - t a 

Iah SH / 


ii ~ *» E x y4 - t a A%g 
t ( A + i 8 ) + 2t ( A + 16) 


(81) 


is continuous at t = 0, at A = 0» & n< l at" ' 9 ™ 0* then 


^(H,0) • / E 1 (H,0,t)f 

-l 


i(A(H) + 10)(1 - t S )]-p4L (82) 

. 2 Vl-t a 


where f is an arbitrary, twice differentiable function of 
one variable will be a solution of 




(83) 


The proof of the above theorem follows step-by-step the 
proof of theorem (63) of reference 3. 

Denote by E 3 (H,0,t) a solution of ^ ( A, -0 , t ) = 0 

and obta.in the following representation for solutions of (83) 
in terms of two arbitrary, twice differentiable functions 
f x , f 3 of one variable. 


>Kh, 0) 



E i ( H , 0,t )fi 
1 


i(A(H) + i 6) (l - t 3 ) 
<2 


— ( A(H ) - 1 6) ( 1 - t*) 
2 


dt 


Jl - t ; 


(84) 


Eor M < 1, A(H) = A(M) (see equation (48) of reference 3) 
is real and therefore £ = A - i 6 , X* A + 10 represent 
conjugate complex variables. For M > 1 , A(H) » ^>(M) be-, 
comes purely imaginary and therefore £ = A + i 8 * i fir + BJ 

= l£ and £ = A_i 0 =i|^A_ 0 ^ = if| where £ and T| are 

the real variables introduced in (64). It remains merely to 
show that (80) can be written in the form (6) for M < 1 and 
in the form (66) for M > 1. Suppose first that M < 1 not- 
ing that $a = s — 5 . = ~ / / ^ 1 (80) can be written in the 

form 3A A H BE' SH 
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/IT 


n 


t 3 L*( * ^-1)1 + i 

H { a, 2 J} tt —pr (_ ^ a02 


a{ 


2A r‘ 


JT 


A 3 a s Bj. BE-! * 

+ a h A HH + BE i 


Since 


a 2 s x 


BE 




3A~ ba 

ft 2 BE . 


0 (85) 


bh 2 \^ba "h j h aA 3 

equation (85) can be replaced, by 


V ♦ Tit A hs 


z/i - t 3 A 


*T 


H / B*® i iL HH BE 
+ 


B £ 4Ajj St 


BBi \ _ 

art / ” t ; 


2 Ah* 


BE x 


£ (v /T-”t S V + 4Ah j 

= 0 (86) 


HH 


E- 


4Ah 3 / B a E i A HH / BEi BEi 

Jl - t»v dJ9 t A E s V a ^ a T 


B 


«V 


B. 



ing factor, it is seen that (87) is essentially the same as 
equation (75) of reference 3. 


Note that t E-; of reference 3 should be corrected to 

-t~ 3* and that (86; is the conjugate of equation (75) of 

reference 3. i.e., the latrer ma.y be obtained from the former 
by replacing £ by and 1[ by £. 
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M 

/* 

1.00 

0.0000 

1.01 

0.0008 

1.02 

0.0020 

1.03 

0.0039 

1.04 

0.0061 

1.05 

0.0085 

1.06 

0.0110 

1.07 

0.0140 

1.08 

0,0169 

1.09 

0.0200 

1.10 

0.0235 

1.11 

0.0268 

1.12 

0.0304 

1.13 

0.0339 

1.14 

0.0376 

1.15 

0.0415 

1.16 

0.0455 

1.17 

0.0494 

1.18 

0.0535 

1.19 

0.0577 

1.20 

0,0618 

1.21 

0.0663 

1.22 

0.0*707 

1.23 

0.0754 

1.24 

0.0^98 

1.25 

0.0845 

1.26 

0,0887 

1.27 

0.0932 

1.28 

0.0981 

1.29 

0.1027 

1.30 

0.1075 

1.31 

0.1124 

1.32 

0.1172 

1.33 

0.1220 

1.34 

0.1268 

1.35 

0.1319 

1.36 

0.1369 

1,37 

0.1417 

1.38 

0.1467 

1.39 

0.1519 

1.40 

0.1567 

1.41 

0.1619 

1.42 

0,1668 

1.43 

0.1719 


0.1772 


0.1819 


0.1873 


M 

/* 

1.47 

0.1921 

1.48 

0.1974 

1.49 

0.2026 

1.50 

0.2078 

1.60 

0.2590 

1.70 

0.3108 

1.80 

0.3618 

1,90 

0.4116 

2.00 

0.4602 

2,10 

0.5076 

2.20 

0.5535 

2.30 

0.5983 

2.40 

0.6413 

2.50 

0.6827 

2.60 

0.7229 

2.70 

0.7613 

2.80 

0.7983 

2.90 

0.8340 

3.00 

0.8682 

3.10 

0.9013 

3.20 

0.9329 

3.30 

0.9637 

3.40 

0.9930 

3.50 

1.0213 

3.60 

1.0487 

3.70 

1.0748 

3-80 

1.1002 

3.90 

1.1243 

4.00 

1.1481 

4.10 

1.1707 

4.20 

1.1926 

4.30 

1.2136 

4.40 

! 1.2339 

4.50 

1.2537 

4.60 

1.2724 

4.70 

1.2908 

4.80 

1.3088 

4.90 

1.3257 

5.00 

1.3424 

5.10 

1.3585 

5.20 

1.3740 

5.30 

1.3891 

5.40 

1.4038 

5.50 

1.4179 

5.60 

1.4315 

5.70 

1.4451 

5.80 

1.4579 


M 

/* 

5.90 

1.4703 

6.00 

1.4827 

6.10 

1.4944 

6.20 

1.5061 

6.30 

1.5172 

6.40 

1.5279 

6.50 

1.5388 

6.6O 

1.5491 

6.70 

1.5591 

6.80 

1.5689 

6.90 

1.5785 

7.00 

1.5877 

7.10 

1.5966 

7.20 

1.6054 

7.30 

1.6143 

7.40 

1.6225 

7.50 

1.6308 

7,60 

1.6388 

7.70 

1.6466 

7.80 

1.6542 

7.90 

1.6617 

8.00 

1.6688 

8.10 

1.6761 

8.20 

1.6829 

8.30 

1.6898 

8.40 

1.6964 

8.50 

1.7028 

8.60 

1.7093 

8.70 

1.7155 

8.80 

1.7216 

8.90 

1.7276 

9.00 

1.7335 

9.10 

1.7391 

9.20 

1.7447 

9.30 

1.7502 

9.40 

1.7557 

9.50 

1.7608 

9.60 

1.7659 

9.70 

1.7710 

9.80 

1.7760 

9.90 

1.7809 

10.00 

1.7858 
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For the case M > 1, a 'similar procedure yields 


2A h exp(- 



i s t£-/L - t 


— Z ( C1 - 




+ a{t 


[ 


3T|3t t 3TI 

S 2 a 1 * 5. / B 

3 £3t) 


_ a- 




88 ) 


which up to a constant factor coincides with (70) . 


APPENDIX IV 


THE COMPLEX POTENTIAL IS THE HODOGRAPH PLANE FOR 

A JOUKOWSKI PROFILE 


1. In connection with the second method for the determi- 
nation of g^(Z) it is necessary to have an analytic rep- 
resentation for the complex potential (in the hodograph plane) 
of an incompressible flow around various profiles. 

This problem will be treated in the following for a sym- 
metric Joukowski profile. 


2. The function 


z + = T)a +- z* + Tie*, T) > 0 


(89) 


maps the circle |z*[ = a into the circle i z + - T|a( = a(l+T]) 
The transformation 


1 / + 4 . 

z = 2 V + TV 


(90) 


maps lz + - Tlal = a(l + T] ) into a Joukowski profile. There- 
fore 


T)a + (1 + T) ) z* + 


T)a + (l + 11) 


(91) 


maps [z*[ = a into a Joukowski profile. 
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Since the complex potential around |z*t = a is 

w(z*) = -V^z , 'e ia + “ 1 log -2— (92) 

V z >n e la / Stt 


z 

a 


the complex potential W(z) = w[z*(z)] is obtained by sub- 
stituting the function 


z* = ( z - 2T)a ) + s 


2(1 + Tl) 

(which is the inverse to (9l))into (92). 
'( z - 2T]a + B)e ia ' 


= ±( z 3 - 4a 3 ) 


¥(z) =* -V 


2a (1 + Tl) 


2(1 + Tl) 


(93) 


> ia ( z - 2Tla + s) J 


IT log 2T ) a -,. + 8 (94) 

2tt 2(1 + Tl)a 


Denoting by 

q 

the 

gives - 



q -3 ,.-18 = M = 

aw 

dz* 

iz 

dz* 

dz 

- r. 

.Ve ia 

+ V 


ir 1 


Q ia z * s 2tt z*J 


+ s + z 

2(1 + Tl)sJ 


(96) 


The aim of this appendix will be to represent V as a 
function of q. By writing 


g) ia ( z - 2T) a + s ) 
2(1 + T]) 


2a (1 + T]) 


ia 


( z - 2T}a + s ) 


“ r j_ ( z, s) (96) 


z - 2Tja + s / v . 

,2(l + T)) a = r s^ z ' 3 ^ (97) 

it is seen that r x and r a are rational functions of z 
and , s, where z, s, and a are connected by the relation 
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5? 


q = -V 


and 


( 


) Ja 


z ) e 


2(1 + 11 ) B 


+ 7 


2(1 + T} )a 2 ( 8 + z) IT 


S + z 


i CC t \ 3 

e (z-2Tla+8)B 


2tr ( z - 2Tla + s ) b 


(98) 


3 3 . 3 

8 = z - 4a 

Introducing a new variable t, defined by 


(99) 


8 ■ (* + t) a 

gives 

8 “ a (‘ - t) 


and it i s found that r x and r 3 become rational functions 
, of t, which will be denoted by and R 2 . . A formal com- 

putation yields 


~(t - Tl)e lCT + (1 + Tl) 

- (1 + T\) (t - ll)eico_ 


R s (t) 


(t - n) 

(1 + Tl) 


t and 
q. - a 3 (t) 


q 


are connect 
0; R s (t) = 


ed by a relation 

e ta U _ P(1 + Tp 

1 + Tl (t-Tj ) s e ia 



+ — — 1 (loo) 

2tTa (t-T))~ 


which 1 b obtained by replacing in (98) e and 
and a (t-ij;, respectively. E i » S a • R 3 


a by a^t+i^) 
are rational 


functions of' t; 2^, Eg are so-called algebraic functions of q, 

The determination of singular points of these functions 
as well as determination of series development of .Rj. and Er a 

around these points can be achieved using classical methods of 
theory of functions. r ~~ "• - -Ti 


The derivation of the corresponding developments of 
and log E 3 as a function of Z = log q does not involve 
additional essential difficulties. 
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Table 4 

The mines of S^and T 
S operation) 


0 = 0=0 


0=0.1 


0=0.3 


S w/ 


0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


T (0) 

s<°> 

x (o) 

s <°> 

T (0) 

s<°> 

T (°) 

s<°> 

0.0000 

-0.0993 

0.0049 

-0.1950 

0.0188 

-0.2841 

0.0394 

-0.3654 

-0,1058 

-0.1113 

-0.0983 

-0.2165 

-0.0777 

-0.3116 

-0.0480 

-0,3955 

-0.2270 

-0.1297 

-0.2152 

-O.2486 

-0.1833 

-0.3509 

-O.I396 

-0.4363 

-0.3735 

-0.1599 

-0.3561 

-0.2982 

-0.3005 

-0.4074 

-0.2344 

-0.4943 

-0.5649 

-0.2147 

-0.5240 

-0.3794 

-0.4303 

-0.4895 

-0.3275 

-O.5643 

-0.8524 

-0.3312 

-0.7491 

-0.5166 

-0,5637 

-0.6062 

-0.4069 

-O.6564 

-1-4440 

-0.6377 

-1.0263 

-0.7424 

-0.6523 

-0.7566 

-0.4445 

-0.7619 

0.0000 

-1.3012 

-0.9458 

-1.0117 

-0.5906 

-0.9102 

-0.4139 

-0.8638 

0.0000 

-1.4350 

-0.3981 

-1.1560 

-0.3995 

-1.0166 

-0.3280 

-0.9449 

0.0000 

-1.2987 

-0.1737 , 

-1.1671 

-0.2397 

-1.0628 

-0.2345 

-0.9943 

0.0000 

-1.2013 

-0.0914 

-1.1355 

-0.1456 

-1.0717 

-0.1617 

-1.0186 


9=0,4 


t (0) 


0=0.5 


g (0) T (0) 


— 0 
-0 

-0.5075 -0.0462 
-0.5573 -0.1086 
-0.6188 -0.1624 
-0.6911 -0.2020 

- 0.2219 
-0.8439 -0.2196 
-0.9198 -0.2015 
-0.9517 -O.1676 
-0.9803 -0.1349 


9=0.6 


g(0) iji(0) 


-0.6099 

-0.6619 
-0.7192 

-0.7791 
-0.8366 

-0.8866 -0.1490 
-0.9260 -0.1323 
-0.9546 -0.1120 


mmrssm 

0=0.8 

0=0.9 

9=1.0 

s (o) 

T (o) 

s<°> 

T (0) 

s<0> 

T (0) 


T (o) 

-0.5643 

0.1383 

' -0.6184 

0.1590 


0.1765 

-0.7136 

0.1910 

-0.5892 

0.0832 

-0.6399 

0.1090 

-0.6859 

0-1308 

-0.7281 

0.1483 

-0.6259 

0.0322 

-0.6660 

0.0629 

-0.7076 

0.0889 

-0.7457 

0.1102 

-O.6565 

-0.0139 

—0*6966 

0.0221 

-0.7328 

0.0517 

-0.7660 

0.0758 

-0.6985 

-0.0525 

-0,7310 

-C.0125 

-0.7683 

0.0201 

-0,7888 

0.0462 

—0.7444 

-0.0818 

-0.7680 

-0.0397 

-0.7908 

-0.0057 

-0.8129 

0,0217 

-0.7916 

-0.1007 ! 

-0.8058 

-0.0587 

-0.8213; 

-0.0249 

-0.8376 

0.0025 

-0.8371 

-0.1090 

-0.8424 

-0.0697 

-0.8510; 

-0.0377 

-0.8619 

-0.0115 

—0.8777 

-0.1080 

-0.8759 

-0.0736 

-0.8788 

-0.0449 

-0.8848 

-0.0209 

-0.9110 

-0.1003 

-0.9048 

-0.0720 

-0.9033 

-0.0474 

-0.9153 

-0.0265 

-0.9381 

-0.0887 

-0.9286 

-0.0666 

-0.9243 

-O.O464 

-0.9237 

-0*0287 
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Table 5 
(S operation)' 

A n,m 



0 

1 

2 

3 

4 

0 

-0.5 

0.62500 

0.119797 

0.024741 

0 . 003871 

H 

Q.g5 

0.31250 

0.1796955 

0.0618525 

0.0135485 

! 

B 

-1 

0.416667 

0.047919 

0 . 007069 

0.000860 

! 

2 1 

-0.6666 

0.166667 

0.013691 

0.001571 

0.000156 

V 

3 

-0.26667 ! 

0 . 047619 

0.003042 

0.000286 

0.000024 

4 

-0.076190 

0.010582 

0.000553 

0.000044 

0.000003 


5 





Table 6 
(S operation) 

^n,m 


x 

0 

1 

2 

3 

4 

i 5 

i . 

1 

0.5708 

0 

0 

0 

0 

. 

0 - 

2 

-0.0817 

+0.5708 

0 

0 

0 

0 

3 

+0.0124 

-0.0817 

+0.2854 

0 

1 

o i 

i 

0 

4 

-0.0016 

+0.0124 

-0.0409 

+0.0951 

0 

0 








































































,-3/2 


1000,00000 

I'm: rwwt 

■»♦*/• WWW 

37.03704 

15.62500 

8.00000 

4.62963 

2.91545 

1.95313 

1.37174 

1.00000 


0.75131 

0.57870 

0.45517 

0.36443 

0.29630 


-1 


100.00000 

nc /ww* 

A./IVWW 

11- 11111 

6.25000 

4.00000 

,2,77778 

2.04082 

1.56250 

1.23457 

1.00000 

0.82644 

0.69444 

0.59172 

0.51020 

O.44444 


Table 7; The Values of p* (0 operation) 


,-V2 


10.00000 

£ n/v w» 

/•VAAAAI 

3.33333 

2.50000 

2.00000 

1.66667 

1.42857 

1.25000 

1.11111 

1.00000 


0.90909 

0.83333 

0.76923 

0.71429 

0.66667 


1. 

*1 

x« 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

1. 

1 . 

1. 

1. 

1 . 

1. 


1/2 


0.1 

A O 

0.3 

0.4 

0.5 

0,6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 


0.01 
A n j 
l/.Uf 

0.09 

0.16 

0.25 

0.36 

0.49 

0.64 

0.81 

1.00 


1.21 

1.44 

1.69 

1.96 

2.25 




0.001 

A AAil 

VI. wo 

0.027 

0.064 

0.125 

0.216 

0.343 

0.512 

0.729 

1.000 

1.331 

1.728 

2.197 

2.744 

3.375 


.2 


0.0001 

G.GGlo 

0.0081 

0.0256 

0.0625 

0.1296 

0.2401 

0.4096 

0.6561 

1.0000 

1.4641 

2.0736 

2.8561 

3.8416 

5.0625 


£ 


>/2 


0.00001 

G.OGC32 

0.00243 

0*01024 

0.03125 

0.07776 

0.16807 

0.32768 

0.59049 

1.00000 

1.61051 

2.48832 

3.71293 

5-37824 

7.59375 


0.000001 

O.GOOG64 

0.000729 

0.004096 

0.015625 

O.O46656 

0.117649 

0.262144 

0.5314a 

1.000000 

1.771561 

2.985984 

4.826809 

7.529536 

II.390625 


->7/2 


0. (6)*1 
0.000013 
0.000219 
0.001638 
0.007813 

0.027994 

0.082354 

0.209715 

0.478297 

1.000000 

1.948717 

3.583181 

6.274852 

10.541350 

17.085938 


f 5 

pll/2 

f 6 

/3/2 

p 7 

^5/2 

0 -( 9 )l 

0.{10)1 

0 .(U )1 

0.(12)1 

0,03)1 

0. (14)1 

0.(6)l 

0.(7)205 

0.(8)410 

0.(9)819 

0.(9)164 

0.(10)328 

nnnnnnfi 

n.nnnm? 

n_ 

n 

n f rA / , 7ft 


0.000105 

0.000042 

0.000017 

0.000007 

0.000003 

0.000001 

0.000977 

0.000488 

0.000244 

0.000122 

0.000061 

0.000031 

0.006047 

0.003628 

0.002177 

0.001306 

0.000784 

0.000470 

0.028248 

0.019773 

0.0138a 

0.009689 

0.006782 

0.004748 

0.107374 

0.085899 

0.068719 

0.054976 

0.043980 

0.035184 

0.348678 

0.313811 

0.282430 

0.254187 

0.228768 

0.205891 

1.000000 

1.000000 

1,000000 

1.000000 

1.000000 

1.000000 

2.593742 

2.853117 

3.138428 

3.452271 

3.797498 

4-177248 

6.191736 

7.430084 

8.916000 

1C.699321 

12.839185 

15-407022 

13-785849 

17.921604 

23.298085 

30.287511 

39.373764 

51.185893 

28.925465 

40.495652 

56.693912 

79.371477 

in.120069 

155-563096 

57-665039 

86,498086 

129.747129 


291.931040 

437.896560 


0.(7)1 

0.000003 

0.000066| 

0.000655 

0.003906 

0.016796 

0.057648 

0.167772 

0.430467 

1.000000 

2.143589 

4.299817 

8.157307 

14.757891 

25.628906 


,9/2 


0 .( 8)1 

0.000001 

0.000020 

0.000262 

0.001953 

0.010078 

0.040354 

0.134218 

0.387420 

1.000000 

2.357948 

5-159780 

10.604499 

20.660X147 

38443359 


o> 

to 


(*) Rote ; The 
ember in pa- 
rentheses in- 
dicates the 
number of zeros 
following the 
decimal point. 
Thus, 0.(7)205 = 
0.0000000205. 
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Table 8: The Values 


<? , 

HU 

HH 

■VI 

m 

1 

oo* 

.000 000 

1.000 000 

.000 000 

1.000 000 

30° 

.258 819 

.965 926 

.500 000 

.866 025 

60° 

.500 000 

.866 025 

.866 025 

.500 000 

90° 

.707 107 

.707 107 

1.000 000 

.000 000 

120° 

.866 025 

.500 000 

.866 025 

-.500 000 

150° 

.965 926 

.258 819 

.500 000 

-.866 025 

180° 

1.000 000 

.000 000 

.000 000 

-1.000 000 

210° 

.965 926 

-.258 819 

-.500 000 

-.866 025 

240° 

.866 025 

-.500 000 

-.866 025 

-.500 000 

•270° 

.707 107 

-.707 107 

-1.000 000 

,000 000 

300° 

.500 000 

-.866 025 

-.866 025 

.500 000 

330° 

.258 819 

-.965 926 

-.500 000 

.866 025 

360° 

.000 000 

-1.000 000 

.000 000 

1.000 000 




% 

cos 5P/2 

sin 3 9 

C08 3 <p 

sin 7f/2 

0° 

1.000 000 

.000 000 

1.000 000 

-000 000 

30° 

.258 819 

1.000 000 

.000 000 

.96$ 926 

60° 

-.866 025 

.000 000 

-1.000 000 

-.500 000 

90 d 

-.707 107 

-1.000 000 

.000 000 

-.707 107 

120° 

.500 000 

.000 000 

1.000 000 

.866 025 

150° 

.965 926 

1,000 000 

.000 000 

.258 819 

IPS 

.000 000 

.000 000 

-1.000 000 

-1.000 000 

210° 

-.965 926 

-1.000 000 

.000 000 

.258 819 

240° 

-.500 000 

.000 000 

1.000 000 

,866 025 

270° 

.707 107 

1.000 000 

.000 000 

-,.707 107 

300° 

.866 025 

.000 000 

-1.000 000 

-.500 000 

330° 

-.258 819 

-1.000 000 

.000 000 

.965 926 

360° 

-1.000 000 

.000 000 

1.000 000 

.000 000 - 


>f cos and sin f(p (0 operation) 


s* 

5 


RH 

■n 

sin 2<p 

■n 

■h 

.000 000 

1,000 000 

.000 000 

1.000 000 

.000 000 

.707 107 

.707 107 

.866 025 

.500 000 

.965 926 

1.000 000 

.000 000 

.866 025 

-.500 000 

.500 000 

.707 107 

-.707 107 

.000 000 

-1.000 000 

-.707 107 

.000 000 

-1.000 000 

-.866 025 

-.500 000 

-.866 025 

-.707 107 

-.707 107 

-.866 025 

.500 000 

.258 819 

-1.000 000 

,000 000 

.000 000 

1.000 000 

1.000 000 

-.707 107 

.707 107 

,866 025 

.500 000 

.258 819 

.000 000 

1.000 000 

.866 025 

-.500.000 

-.866 025 

.707 107 

.707 107 

.000 000 

-1.000 000 

-.707 107 

1.000 000 

.000 000 

-.866 025 

-.500 000 

.500 000 

.707 107 

-.707 107 

-.866 025 

.500 000 

.965 926 

.000 000 

-1.000 000 

.000 000 

1.000 000 

.000 000 


cos 7f/2 

sin 

49 

cos 

49 

sin 9f/2 

COO 9^/2 

1,000 

000 


m 

1.000 

000 

.000 

000 

1.000 

000 

-«2 58 

819 


025 

-.500 

000 

.707 

107 

-.707 

107 

-.866 

025 


025 

-.500 

000 

- 1.000 

000 

.000 

000 

.707 

107 

.000 

000 

1.000 

000 

.707 

107 

.707 

107 

.500 

000 

.866 

025 

-.500 

000 

.000 

000 

- 1.000 

000 

-965 

926 

-.866 

025 

-.500 

000 

-.707 

107 

.707 

107 

.goo 

000 

.000 

000 

1.000 

000 

1.000 

000 

.000 

000 

.965 

926 

.866 

025 

-.500 

000 

-.707 

107 

-.707 

107 

-.500 

000 

-.866 

025 

-.500 

000 

.000 

000 

1.000 

000 

-.707 

107 

.000 

000 

1.000 

000 

.707 

107 

-.707 

107 

.866 

025 

.866 

025 

-.500 

000 

- 1.000 

000 

.000 

000 

.258 

819 

-.866 

025 

-.500 

000 

.707 

107 

.707 

107 

- 1.000 

000 

-000 

000 

1.000 

000 

.000 

000 

i.opo 

000 


•* 


05 

03 
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Table 9 (S operation) ^ ACA No * 
Computation of the stream function (in the logarithmic plane) of 


1018 


■ 

6 

s<o> 


g (2) 

s (3) 

s (4) 

-£S_ 

»A 

7 ) A 

.02 

.05 

-.0496 


-.0002 

^.0004 

.0001 

.0538 



.2 

-.1945 


.0013 

.0000 

.0001 

.1862 



.4 

-.3627 


.0105 

-.0001 

-.0001 

.2653 



.6 

-.4970 

-.0172 

.0342 

-.0003 

-.0009 

.2511 



.8 

-.6223 

-.1035 

.1764 

.0003 

-.0034 

.2125 

.2006 

.06 

.08 

-.0840 

.0047 

,0000 

.0000 


.0984 

.0999 


.30 

-.2961 

.0142 

.0033 

-.0005 


.2792 

.2835 


.34 

-.3147 

.0112 

.0058 

-.0001 


.2817 

.2975 


.40 

-.3584 

.00*34 

.0091 

-.0004 


.2906 

.3085 


.50 

-.4460 

.0096 

.0200 

-.0011 


.2996 

.3056 


.70 

-.5640 

.0054 

.0492 

-.0043 

-.0020 

.2575 

.2541 


.90 

-.6577 

*.0322 

.1083 

-.0020 

.0109 

.1927 

.1832 

.10 

.15 

-.1621 

.0148 

-.0001 

.0000 

*.0000 

.2050 

.2082 


.35 

-.3490 

.0275 

.0059 

-.0066 

-.0000 

.3380 

.3424 


.50 

-.4593 

.0282 

.0194 

-.0018 

-.0003 

.3362 

.3406 


.60 

-.5218 

.0246 

.0338 

-.0030 

-.0008 

.3109 

.3134 


.70 

-.5770 

.0178 

.0539 

-.0044 

-.0016 

.2529 

.2773 


1.00 

-.7085 

-.0179 

.1507 

-.0094 

-.0083 

.1993 

.1608 

.20 

.22 

-.2672 

.0458 

-.0024 

.0000 

.0000 

>4119 

.4150 


.40 

-.4311 

.0706 

.0039 

-.0017 

.0001 

.4701 

.4731 


.70 

-.6145 

.0802 

.0623 

-.0057 

-.0010 

.3224 

.3549 


.80 

-.6556 

.0722 

.0730 

-.0093 

-.0017 

.2907 

.2842 


1.00 

-.7324 

.0553 

.1495 

-.0244 

-.0059 

.2092 

• .1908 


1.10 

-.7664 

.0442 

'.1971 

-.0301 

-.0101 

.1781 

.1508 

.30 

.30 

-.4454 

.1193 

-.0143 

.0005 

.0001 

.9389 

.6771 


.75- 

-.6934 

.1791 

.0422 

-.0197 

.0015 

.3896 

.3593 


.85 

-.7274 

.1791 

.0726 

-.0288 

.0014 

.3244 

.2959 


.95 

-.7602 

.1757 

.1102 

-;0396 

.0006 

.2648 

.2411 


1.10 

-.7158 

,1666 

.1814 

-.0589 

-.0025 

.2132 

.1722 


1.20 

-.8317 

.1578 

.2388 

-.0740 

-.0065 

.1411 

.1334 

.40 

.35 

-.5895 

.1873 

-.0326 

.0019 

.0001 

1.0839 

.9059 


.60 

-.6962 

.2439 

-.0128 

-.0085 

.0021 

.5829 

.5458 


.85 

-.7672 

.2587 

.0493 

-.0351 

.0049 

•3458 

.3.240 


.93 

-.7883 

.2588 

.0789 

-.0469 

.0057 

.2980 

.2737 


1.05 

-.8186 

.2553 

.1333 

-.0673 

.0054 

.2416 

.2109 


1.30 

-.8764 

.2197 

.2881 

-.0926 

-.0074 

.1639 

.1122 

.60 

.40 

-.7621 

.3008 

-.0653 

.0074 

-.0001 

1.0761 



.45 

-.7654 

.3149 

-.0652 

.0055 

.0005 

.9092 

.9117 


.60 

-.7808 

.3460 

-.0532 

-.0042 

.0036 

.599C 

.5959 


.75 

-.8028 

.3628 

-.0240 

-.0227 

.0072 

.4250 

.4188 


.95 

-.8378 

.3714 

.0434 

-.0595 

• 0133 

.2869 

.2734 


1.00 

-.8479 

.3718 

,0653 

-.0709 

.0147 

.2623 



1.15 

-.8783 

.3694 

.1440 

-.1108 

.0180 

.2045 

.1785 


(continued on next page) 
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Table 9 (Continued) 


-A 

9 

HO 

H 

811 

HI 

1 

££££ 

>8 
b X 

.80 

.50 

-.9096 

.4958 

-.1450 

.0233 

-.0012 

.5409 

.0425 


.70 

-.8844 

.5258 

-.1238 

-.0005 

.0076 

.3759 

.3753 


.90 

-.8924 

.5407 

-.0676 

-.0453 

.0212 

.2666 

.2625 


1.05 

-.9090 

.5448 

-.0020 

-.0928 

.0334 

.2089 

.1999 


1.10 

-.9166 

.5449 

.0236 

-.1129 

.0396 

.1916 

.1819 


1.20 

-.9323 

.5445 

.0832 

-.1526 

.0457 

.1648 

.1496 


1.40 

-.9721 

.5401 

.2281 

-.2515 

.0592 

.1242 

.0959 

1.00 

• 54 

-.9729 

.6908 

-.2628 

.0602 

-.0075 

.2251 

.2292 


.75 

-.9428 

.7120 

-.2373 

.0246 

.0093 

.2156 

.2211 


1.00 

-.9453 

.7247 

-.1534 

-.0589 

.0426 

.1612 

.1685 


1.10 

-.9538 

.7270 

-.1036 

-.1045 

.0586 

.1319 

.1459 


1.20 

-.9664 

.7277 

-.0446 

-.1578 

.0756 

.1191 

.1246 


1.35 

-.9903 

.7262 

.0611 

-.2522 

.1020 

.0927 

.0949 

1.20 

.58 

-.9925 

.8907 

-.4198 

.1231 

-.0231 

.0816 

.0892 


. .80 

-.9726 

.9049 

-.3869 

.0722 

.0070 

.1082 

.1230 


1.05 

-.9748 

.9131 

-.2937 

-.0397 

.0623 

.0869 

.1108 


1.15 

-.9830 ■ 

.9152 

-.2405 

-.1001 

.0895 

.0725 

.1000 


1.22 

-.9899 

.9157 

-.1945 

-.1495 

.1098 

.0625 ' 

.0919 


1.30 

-1.0008 

.9149 

-.1411 

— 2090 

.1346 

.0508 

.0822 


1.40 

-1.0171 

.9159 

-.0653 

-.2924 

.1662 

.0357 

.0699 


(continued on next page) 


Note; Columns 1 - 6^, i.e. were computed "by means of 

series as described in method II j column. 7 namely was 

computed directly from the formula for the stream function, see refer- 
ence 3 • 













Table 9 


g j(0) T (l) j(2) t ( 3) »p(A) 


.05 

mmti 

■KlSfri 

MKOiM 

B$if| 

-.0014 


,0003 

-,0198 

,0004 

nrwrf 

.40 

.0476 

-.0773 

.0006 

.0012 

.0002 

.60 

.1005 

-.1658 

-.0010 

.0054 

.0002 

.80 

.1503 

-.2627 

-.0332 

.0156 

-.0000(01) 

.08 

-.0576 

-.0015 

.0002 

.00001 

.0000(4) 

.30 

-.0129 

-.0443 

.0018 

.00003 

.0002 

.34 

.0006 

-.0533 

-.0010 

.0012 

.0000(5) 

.40 

.0195 

-.0719 

-.0034 

.0019 

.0000(2) 

.50 

.0468 

-.1204 

.0050 

.0027 

-.0001 

.70 

.0917 

-.2160 

.0080 

.0086 

-.0008 

.90 

.1390 

-.3407 

.0020 

.0236 

-.0014 

1.10 

.1565 

-.4689 

-.0008 

.0478 

-.0070 

.35 

-.3338 

.0123 

.0188 

-.0048 

.0006 

.60 

-.1228 

-.1559 

.0725 

-.0112 

.0001 

.85 

-.0100 

-.3397 

.1366 

-.0074 

-.0053 

.93 

.0152 

-.4023 

.1582 

-.0023 

-.0086 

1.05 

.0476 

-.4978 

.1879 

.0106 

-.0156 

1.30 

.1650 

-.7081 

.2173 

.2080 

-.0584 

.15 

-.0893 

-.0073 

.0010 

-.00002 

.0001 

.35 

-.0317 

-.0598 

.0054 

.0005 

.0001 

.50 

.0198 

-.1218 

.0098 

.0023 

-.0001 

.60 

.0519 

-.1717 

.0119 

.0048 

-.0003 

.70 

.0816 

-.2282 

.0145 

.0091 

-.0008 

1.00 

.1488 

-.4270 

.0147 

.0376 

-.0039 

.22 

-.1753 

-.0091 

.0035 

-.0014 

.0019 

.40 

-.0923 

-.0730 

.0141 

-.0006 

.0000(4) 

.70 

.0604 

-.2573 

.0326 

.0095 

-.0005 

.80 

.0626 

-.2984 

.0453 

.0127 

-.0030 

1.0Q 

.1123 

-.4402 

.0588 

.0347 

-.0075 

1.10 

.1316 

-.5165 

.0631 

.0513 

-.0108 



8T0T -o/i KI TOVit 











tl- ' - * - f I ^ > — J u \ v »rn — 

described in method XI> columns 7-11, l.e., T^°>, ..., yf were ipomputed 

directly from the formula for the stream functions, see reference 3. 


Table 9 


ST -X 0 
•* _ 


t (2) t (3) t (4) 


.30 

.30 

-.2856 

-.0012 

.0111 

-.0019 

-.0108 


.75 

-.0165 

-2681 

.0821 

-.0014 

—0069 


.85 

+.0208 

-3396 

.1008 

+.0053 

—0056 


.95 

+.05&1 

-4142 

.1183 

+.OH44 

-.0092 


1.10 

+.0855 

-5320 

.1438 

+.0358 

-.0169 


1.20 

+.0985 

-.6144 

.1600 

+.0566 

-.0240 

.60 

.40 

-3132 

+.OI68 

.0291 

-.0102 

+.0019 


•45 

-.2639 

—0214 

.0447 

-.0135 

+.0022 


.60 

-1533 

-1369 

.0943 

+.0224 

+.0015 


.75 

—0777 

-2552 

.1475 

+.0286 

+.0004 


.95' 

-0076 

-4180 

.2212 

+.0270 

-.0075 


1.00 

+.0060 

-4597 

.2398 

+.0244 

-.0107 


1.15 

+.0406 

-.5872 

.2955 

+.0088 

-.0240 

.80 

.50 

-.1986 

-0163 

.0679 

-0297 

+.0070 


.70 

-.1078 

-.1939 

.1704 

-.0571 

+.0097 


.90 

—.0446 

-.3692 

.2775 

-.0767 

+.0058 


1.05 

—0094 

-.5014 

.3587 

-.0817 

—0041 


1.10 

+.0007 

-.5463 

.3864 

-.0814 

—0090 


1.20 

+.0183 

-.6358 

.4408 

-.0757 

-.0216 


1.40 

+.0465 

-.8183 

.5505 

-.0447 

-.0596 

.00 

.54 

—0619 

+.3863 

-.5576 

-0538 

+.0168 


.75 

+.0923 

+.3747 

-.7101 

—1071 

+.0263 


1.00 

+.2105 

+.3884 

-9227 

-.1568 

+.0235 


1.10 

+.2556 

+.4009 

-1.0174 

-.1697 

+.0157 


1.20 

+.2995 

+.4164 

-1.1175 

-.1767 

+.0033 


1.35 

+.3640 

+.4450 

-1.2776 

-1752 

-0255 

..20 

.58 

+.0785 

+.4983 

-7305 

—0903 

+.0338 


.80 

+.1662 

+.5113 

—8962 

-1796 

+.0566 


1.05 

+.2708 

+.5484 

-1.1099 

-.2664 

+.0621 


1.15 

+.3122 

+.5696 

-1.2038 

-2934 

+.0556 


1.22 

+.3412 

+.5857 

-1.2717 

-.3110 

+.0475 


1.30 

+.3733 

+.606I 

-1.3535 

-3215 

+.0335 

_ 

1.40 

+,4133 

+.6350 

-1.4595 

-.3322 

* 

+.0102 



-.1260 

-.0780 

-.0194 

-.0019 

+.0007 

+.0171 

-.0757 

-.0511 

-.0213 

-.0107 

-.0040 



-.0159 

-.1939 

-.3689 

-.5020 

-.5467 

-.6362 

-.8187 

-.0218 

-.2213 

-4525 

-5451 

-.6376 

—7782 

-0414 

-.2558 

-.4945 

-.5896 

-.6563 

-7325 

—8283 


HACA 
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Table 10 

The values of . A * 

AA A A A0 AO 



A =0 

g»K£ 

SH 

A = 

= 0.2 

>- 0.3 

A = 0.4 

MKSR91 

o- 

& s (°) 


AS(0) 

at (°) 

AS(0) 

at (°) 

^ s ( O ) 

at(°) 

4S(°> 

At (°) 

as(O) 


A A 

A A 

Aa 

A A 

AA 

4A 

4A 

A* 



' 

0 

.0 

d . 6o 

0.7 

0.00 

1.05 

0.00 

1.40 

0.00 

1.75 

0.00 

2.1 

0.00 


0 

.1 

0.13 

0.8 

0.15 

1.05 

0.17 

1.30 

0.46 

1.5 

0.75 

1.7 

1.75 

3.0 

0 

.2 

0.2 

0.9 

0.33 

1.07 

0.46 

1.24 

0.6 

1.4 

1.2 

1.3 

1.8 

1.2 

0 

.3 

0.2 

0.9 

0.33 

0.9 

0.46 

0.9 

0.6 

0.9 

1.0 

0.8 

1.4 

0.7 

0 

.4 

0.3 

0.65 

0.4 

0.70 

0.5 

0.75 

0.6 

0.8 

0.85 

0.625 

1.1 

0.45 

0 

.5 

0.33 

0.7 

0.42 

0.65 

0.51 

0.59 

0.6 

0.53 

0.675 

0,415 

0.75 

0.3 

0 

.6 

0,28 

0.6 

0.35 

0.56 

0.42 

0.53 

0.5 

o.s 

0.55 

0.365 

0.6 

0.25 

0 

.7 

0.23 

0.55 

0.29 

0.49 

0.35 

0.43 

0.4 

0.37 

0.43 

0.31 

0.46 

0.248 

0 

.8 

0.18 

0.5 

0.22 

0.46 

0.26 

0.42 

0.3 

0.38 

0.327 

0.31 

0.354 

0 . 24 1 

E 


0.14 

0.44 

0.18 

0.40 

0.22 

0.36 

0.26 

0.33 

0.273 

0.28 

0.286 

0.223 

B 

E 

0.1 

0.42 

0.14 

0.385 

0.18 

0.35 

0.21 

0.315 

0.25 

0.28 

0.246 

0.22 



* Note t^LtS§i£l = jgl°l ,4£ ( -°2 = - AS ( 0 ). . 

a> £© &© 

(S operation) 
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69 


^(X, ©) in the vicii 

(S oi 

iity of the curve i|r = 0. 
peration) 

A 

9 

i 


A 

9 

V 

i 

• 

o 

JO 

.10 

-.014 

t .60 

.90 

-.018 


.20 ■ 

.006 


.95 

-.003 


• 30 

.038 


1.00 

.009 





1.05 

.020 

- ,06 

.20 

-.036 





.30 

.003 

‘ - .80 

.90 

-.040 


.40 

.038 


1.05 

-.004 





1.10 

.007 

- .10 

.30 

-.043 


1.20 

.025 


.40 

-.001 


1.40 

.055 


.45 

.017 





.50 

.036 

-1.00 

1.00 

-.026 





1.05 

-.018 

- .20 

.40 

-.088 


1.10 

-.010 


.50 

-.038 


1.15 

-.003 


.55 

-.014 


1.20 

.004 


.60 

.007 



(► 


.65 

* .026 

-1.20 

1.10 

-.014 





1.20 

-.004 

i 

• 

u> 

o 

.65 

-.024 


1.30. 

.005 


.70 

.001 


1.40 

.013 


.75 

.018 





.80 

.032 

-1.90 

1.3 

-.0025 





1.4 

-.0002 

- .40 

.70 

-.033- 


1.45 

.0010 


.80 

-.001 


1.5 

.0020 


.85 

.011 





.95 

.032 





Table 12 


The values of a^y, Yq along the ourve 1jr = 0. 


(S operation) 


A 

.A 




1.024 

1.085 . 

.256 


1.055 

.990 

.377 


1.002 

.845 

.400 

- .20 

.852 

.580 

.431 

- .30 

.715 

.410 

.392 

- .40 

.528 

.260 

.286 

- .50 

.483 

.208 

.274 

- .60 

.418 

.159 

.269 

- .80 

.366 

p .110 

.229 

-1.00 

.254 

.061 

.134 

-1.20 

.166 

.031 

.090 

-1.90 

.067 

,006 

.023 
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Table 13 
(S operation) 

© as a function of v/a Q for one quadrant of the curve *P = 0 


—I 

v/a 0 

M 

T 

A 

© 

.707 

.74523 

.66681 

-0.12091 

.177 

.669 

.70111 

.71305 

-0.15905 

.307 

.633 

.66000 

.75127 

-0.19962 

.403 

.557 

.57513 

.81806 

-0.30026 

.582 

.494 

.50652 

.86223 

-0.40123 

.698 

.440 

.44877 

.89365 

-0.50296 

.805 

.395 

.40131 

.91594 

-0.60065 

.885 

.355 

.35956 

.93312 

-0.69930 

.962 

.288 

.29042 

.95690 

-0.89731 

1.067 

.234 

.23529 

.97193 

■ -1.09784 

1.168 

.183 

.18362 

.98300 

-1.33821 

1.245 

.089 

.08907 

.99603 

-2.05277 

1.406 
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Table 14 
(S operation) 

fit — 8» 

The values of —=■ x, — y along the curve Y = 0, for one quadrant cf the 
ns Po 

curve. The curve is symmetric with respect to both the x and y axes. 


v /*o 

M 

T 

A 

>lrl 

r O 

'i 1 1 1 ' n 

.10 

0.10010 

•99498 

-1.93676 

.027 

.002 

.20 

0.20081 

.97963 

-1.25082 

.055 

.009 

• 30 

0.30274 

•95307 

-0.85821 

.080 

.020 

.40 

0.40656 

.91362 

-0.58907 

.103 

.034 

.50 

0.51299 

.85839 

-0.39083 

.123 

.052 

.55 

0.56743 

.82342 

-0.31062 

.131 

.063 

.60 

0.62284 

.78235 

-0.24069 

.138 

.077 

.65 

0.67933 

.73383 

-0.17949 

.146 

.093 - 

.70 

0.73704 

.67585 

-0.12761 

.152 

.113 

.725 

0.76640 

.64236 

-0.10237 

.156 

.135 


Table 15 - 

(S operation) 

The values of ~ ilZ > along the curve / Y r = 0 for one quadrant of the 

t o dv Po dv 


curve. 


M 

v/a Q 

T 

■* 1 1 - 1 

A 

~ a o dv 
/o dv 

“ a o dx 

Po dv 

.766 

.725 

.64284 

-0.10477 



.745 

.707 

.66707 

-0.12111 

.99 

5.54 

.701 

.669 

. .71316 

-0.15921 

1.29 

4.07 

,660 

.633 

.75127 

-0.19961 

1.48 

3.49 

.575 

.557 

.81816 

-0.30044 

1.69 

2.56 

.507 

.494 

.86195 

-0.40029 

1.82 

2.17 

.401 

.395 

.91608 

-0.60128 

1.95 

1.59 

.360 

.355 

.93295 

-0.69820 

2.19 

1.53 

.235 

.234 

.97200 

-1.09880 

2.65 

1.13 

.184 

.183 

.98292 

-1.33648 

2.75 

.94 

.089 

.089 

.99603 

-2.05277 

- 2.80 

.45 
























Figure 5.- The hodograph of a flow around the profile in fig. 4 
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| Figure 10,- Tbfl value* of B^(*> .3), T^{A, ,3), 


Figure. 11.- Ihi Ttlue* at ■(“id, .1), ?W(x, .1). 
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figure IS.- the values of ^°5(-.3,0), t^ o) (--3.e). 
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figure 14.- The Inga of X*,B) ■ 0 in the phytloal plana. 

{The oontour of the oosrpreaeible flow obtained 

froa the funotion g(2) - 1/2 ti-Se 2 ) 1 ^ 3 + (l-Se 2 )' 1 ^ 3 ], 

® \ t iO) . 



figure 1*.- The valuea of y* , v 5 . 

po dy ' t 0 dT 


> 04 






SAGA TH SO. 1018 



Figure 
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Fige. 18,1? 



.861 .821 U 

16.- The function* T and T 10 . 



coefficients -o. n ; bn. 





